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PREFACE 


This book was undertaken at the suggestion of my lamented 
colleague Professor Benjamin Osgood Peirce, and with the promise 
of his collaboration. His untimely death deprived me of his invalu- 
able assistance while the second chapter of the work was still 
unfinished, and I have been obliged to complete my task without 
the aid of his remarkably wide and accurate knowledge of Mathe- 
matical Physics. 

The books to which I am most indebted in preparing this treatise 
are Thomson and Tait’s “Treatise on Natural Philosophy,’ Watson 
and Burbury’s “ Generalized Coordinates,” Clerk Maxwell’s ‘ Elec- 
tricity and Magnetism,” E. J. Routh’s “Dynamics of a Rigid 
Body,” A. G. Webster’s “ Dynamics,” and E. B. Wilson’s “ Advanced 
Calculus.” 

For their kindness in reading and criticizing my manuscript I am 
indebted to my friends Professor Arthur Gordon Webster, Professor 
Percy Bridgman, and Professer Harvey Newton Davis. 


W. E. BYERLY 
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GENERALIZED COORDINATES 


CHAPTER I 
INTRODUCTION 


1. Coordinates of a Point. The position of a moving particie 
may be given at any time by giving its rectangular codrdinates 
x, y, 2 referred to a set of rectangular axes fixed in space. It 
may be given equally well by giving the values of any three 
specified functions of 2, y, and zg, if from the values in question 
the corresponding values of 2, y, and z may be obtained uniquely. 
These functions may be used as coérdinates of the point, and 
the values of 2, y, and z expressed explicitly in terms of them 
serve as formulas for transformation from the rectangular sys- 
tem to the new system. Familiar examples are polar coordi- 
nates in a plane, and cylindrical and spherical coordinates in 
space, the formulas for transformation of codrdinates being 


respectively F , 

Z =f COs >, Tia COS O; 
| Seal (1) eee, (2) y=rsin@cosd¢,; (3) 
Lie | 22 | z=rsn@singd 


It is clear that the number of possible systems of codrdinates 
is unlimited. It is also clear that if the point is unrestricted in 
its motion, three codrdinates are required to determine it. If it 
is restricted to moving in a plane, since that plane may be taken 
as one of the rectangular codrdinate planes, two coordinates are 
required. 

The number of independent codrdinates required to fix the 
position of a particle moving under any given conditions is 


called the number of degrees of freedom of the particle, and is 
1 
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equal to the number of independent conditions required to fix 
the point. 

Obviously these coordinates must be numerous enough to fix 
the position without ambiguity and not so numerous as to render 
it impossible to change any one at pleasure without changing 
any of the others and without violating the restrictions of the 
problem. 

2. Dynamics of a Particle. Free Motion. The differential 


equations for the motion of a particle under any forces when 
we use rectangular codrdinates are known to be 


ML= KX, 
mij = ¥. (ys 
Die= Z. 


X, Y, and Z, the components of the actual forces on the particle 
resolved parallel to the fixed rectangular axes, or rather their 
equivalents mz, my, mz, are called the effective forces on the parti- 
cle. They are of course a set of forces mechanically equivalent 
to the actual forces acting on the particle. 

The equations of motion of the particle in terms of any other 
system of codrdinates are easily obtained. 

Let ¢,, 9) Y,. be the codrdinates in question. The appropriate 
formulas for transformation of coordinates express x, y, and z in 
terms of q,, q,, and Is 


r=f,(d» I) q3)> y =f, Ga Is)» 2=f,(4) q 93) 
For the component velocity « we have 


ee, Hee» aes 
Oa Ogu 09, a 
and #, 7, and 2 are explicit functions of ¢,, 9) G4. 9» J,» and dy 
linear and homogeneous in terms of 4,, ,, and 4,. 





* For time derivatives we shall use the Newtonian flucion notation, so that 


2 
we shall write « for Us % for de 
dt dt2 
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We may note in passing that it follows from this fact that 
z, gy’, and 2 are homogeneous quadratic functions of Giada 
and 4,. 








5 hare 
Obviously ae _ oe, es 
is eg, 
and since d or x. ae Ox eae 
: - Seyret - - 
dt eq, ag” bq,69, "4,09, * 
6c Ox Ona ox 
and m : 
oq, og; q, ace 4, €4, Y, ar 04,84, G3) 
d ox Ox 


dts =o S 


Let us find now an expression for the work 6, W done by the 
effective forces when the coordinate gq, is changed by an infini- 
tesimal amount 6g, without changing g, or q,. If 62, dy, and 
dz are the changes thus produced in 2, y, and z, obviously 


5,,.V = m[adx + yoy + 262] 
if expressed in rectangular codrdinates. We need, however, to 
express 5, W in terms of our coordinates q,, g,, and q,. 


OL Oye. OF 
6, V= 6g.. 
Is m E 24, ae 0) bq, + a | 4 








oe Cf. Ge . a Ox 
Sy e oq, ~ dt (2 ie eee dt 29,’ 
but by (2) and (3) 
6c te | don 


ee IN 
eq, 4, dt oq, eq, 


o Boe eh tye 2 B (eaiearta( 2) 
pond a Og, dt 0g,\2/) a9,\2)? 





Cpe “atid i 
VN | Se ee 6 Be 4 
and therefore 6,,W | di 04, al G (4) 
where LT = 4 [Pty +2] 


and is the kinetic energy of the particle. 


4 INTRODUCTION [Arr. 3 


To get our differential equation we have only to write the 
second member of (4) equal to the work done by the actual 
forces when q, is changed by 6q,. 

If we represent the work in question by @,6q,, our equation is 


oh OH OF 
eee ee Oe 5 
dt oq, 4, ®, o> 


and of course we get such an equation for every codrdinate. 

It must be noted that usually equation (5) will contain gq, 
and q, and their time derivatives as well as OF and therefore can- 
not be solved without the aid of the other equations of the set. 

In any concrete problem, 7 must be expressed in terms of q,, 
J. Vy and their time derivatives before we can form the expres- 
sion for the work dene by the effective forces. @,69, ,64,) 
(),0q,, the work done by the actual forces, must be vbtained 
from direct examination of the problem. 


3. (a) As an example let us get the equations in polar codr- 
dinates for motion in a plane. 


Here Z=P7 COS , y=rsin >. 
a ae Vi = vy? = 7 as Os 
ER pO 
and of oars Es +7"). 
— = mr 
or ; 


I ss 
Ca mrp. 

8,.W = m[7 —rdé”] dr = RSr 

if # is the impressed force resolved along the radius vector. 

Es = mr, 
ar 
op 
dia 

5,V =m Fi (7°p) 86 = Brdh 


0. 
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if ® is the impressed force resolved perpendicular to the radius 
vector. 
In more familiar form 


ete ils 
ie | (o ©) 


sue 
Poiudt\. ‘dip 


(6) In cylindrical codrdinates where 


Z=7 COs d, y=rsin d, Z= 8, 


T= 5 [P+ rd? +2]. 


oT : 
sa 
or 

or , 
—=mr 
or %% 
oe — mr, 
op 

Cae : 
=~ = ™mMé. 
CZ 


SW =m[i ~ r¢’] dr = Ror, 
d 25 o 
5,VW =m a (7°) df = Ord, 


5,W = m2zb2z = Z0z ; 


d’r ay |- 
lines R, 
BP ie Iz (of 3 
md (ot). 
rd dt = 
D) 
m giz = Z, 
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(c) In spherical codrdinates where 
Z=r cos U, y=r sin @ cos q, z=rsin @sin 9, 
T= - [7? + 7°68? + 7° sin? 66"). 
os mer, 
a 


= = mr [ee ++ sin? 647}, 
la 


— = mr’é, 
é 
— = mr’ sin 6 cos 6¢%, 


ye = mr? sin? Od. 
ag 
6W=m[r—r (C + sin? 04") ] or = Rér, 


5,V =m I< (7°60) —r* sin 6 cos a4 | 50 = Ordé, 
os = me (7 sin? 06) op = Pr sin 489 ; 
ae (a) +) |= 
or Pag ‘ + sin’é ap ae 


mle (et) na ee (2) |- 
ml Sr FE r’ sin 6 cos 0 7 = @, 


ae S(r sin? 6 =)- ®. 


r sin 6 dt 





4. Dynamics of a Particle. Constrained Motion. If the particle 
is constrained to move on some given surface, any_two inde- 
pendent specified functions of its rectangular codrdinates a, y, z, 
may be taken as its coordinates qg, and g¢,, provided that by the 
equation of the given surface in rectangular codrdinates and 
the equations formed by writing q, and q, equal to their values 
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in terms of 2, y, and z the last-named codrdinates may be 
uniquely obtained as explicit functions of q, and g,. For when 
this is done, the reasoning of Art. 2 will hold good. 

If the particle is constrained to move in a given path, any 
specified function of 2, y, and z may be taken as its codrdinate 
4, provided that by the two rectangular equations of its path 
and the equation formed by writing g, equal to its value in 
terms of 2, y, and z the last-named codrdinates may be uniquely 
obtained as explicit functions of g,. For when this is done, the 
reasoning of Art. 2 will hold good. 


5. (a) For example, let a particle of mass m, constrained to 
move on a smooth horizontal circle of radius a, be given an 
initial velocity V, and let it be resisted by the air with a force 
proportional to the square of its velocity. 

Here we have one degree of freedom. Let us take as our 
codrdinate g, the angle @ which the particle has described about 
the center of its path in the time ¢ 


M of 
Yb = my a’@?, 

oT Mh 

and we have = mad, 





oo 
Our differential equation is 
mebs0 = — ka?6?a8d0, 


which reduces to ee . al? = 0, 
or qe -f- ie ab? = 0. 
dt m 


Separating the variables, 


ge + ug adt = 0. 
62 =m 
a 


Integrating, _ ; 4 . at= C=— a 
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le ma+kVat 
6 mV 
dO mV 

dt ~ ma+kVat’ 





= — log [m+kVt] +6, 


ka 
m kVt 


and the problem of the motion is completely solved. 

(6) If, however, we are interested in A, the pressure of the 
constraining curve, we must proceed somewhat differently. We 
have only to replace the constraint by a force & directed toward 
the center of the path. There are now two degrees of freedom, 
and we shall take @ and the radius vector 7 as our coordinates 
and form two differential equations of motion. 


aa) m ey) 242 
i oy [7 + 76 ]. 


i = mr’é, 
ge = mr 
fo = 
= — mr? 
- 
d 20 3A2 
mM °F (770) 50 = — kr? 80. (1) 
m(r — 70”) or =— Rbr. (2) 
To these we may add 
r= da. 
Whence O+ “e 6? = 0, as before, (3) 


and R = mae”. (4) 


Cuap. I] CONSTRAINED MOTION OF A PARTICLE 9 


(¢) Let us now suppose that the constraining circle is rough. 
Here, since the friction is w (the coefficient of friction) multi- 
plied by the normal pressure, R will be needed, and we must 
replace the constraint by & as before. 

We have now 


mn ‘ (7°6) 80 = — kr°6°30 — pRrd6, 
mr — 6") or = — Rr, 


and f=. 
Whence R = ma6*, as before. 


6+ oa. we 0, 


safe as|eao 


Replacing “ in Art. 5, (a), (1), by + ay 
t 
log E + (= + H) “|. GL) 





we have 6= 


a 
— +p 
a 


EXAMPLES 


1. Obtain the familiar equation as see sin6=0 for the 
simple pendulum. 
2. Find the tension of the string in the simple pendulum. 


2 
Ans. k= ml g cos 6 + (5) | 
dt 
3. Obtain the equations of the spherical pendulum in terms 
of the spherical coordinates @ and ¢. 
Ans. 6 — sin 6 cos 6¢° . sinvd.==(), sin? 06 = C. 


6. (a) The constraint may not be so simple as that imposed 
by compelling the moving particle to remain on a given surface 
or on a given curve. 
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Take, for example, the tractriz problem, when the particle 
moyes on a smooth horizontal plane. 

Let a particle of mass m, attached to a string of length a, rest 
on a smooth horizontal plane. The string lies straight on the 
plane at the start, and then the end not attached to the particle 
is drawn with uniform velocity along a straight line perpendicu- 
lar to the initial position of the strimg and lying in the plane. 
Let us take as our codrdinates 2, the distance traveled by 
that end of the string which is not attached to the particle, 
and @, the angle made by the string with its initial position. 
Let & be the tension of the string and » the velocity with 
which the end of the string is drawn along. Let X, Y, be the 
rectangular codrdinates of the particle, referred to the fixed 
line and to the initial position of the 
string as axes. 


X=2r—<asin 0, 
Y=acos?; 


X=%—acos 06, 





Y=— asin 06. + 


=F(24+¥)= Sle +0 — 2a cos 046). 


oT ; ; 
re m[& — a cos 66], 


a = m[a6 — a cos 6<], 
a = ma sin 626. 


00 
ap . 5 
m oF [z — a cos 00] x = R sin 062. 


™ E (a6 — a cos 6%) — asin oid | 60 = 0, 
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Adding the condition sean, 


and reducing, —ma[cos 06 — sin 667]= RF sin 0. 


ma?é = 0. 
6=0. 
R= mab, 
Integrating, 6=C=*. 
a 
ee mn’ 
pana 


The particle revolves with uniform angular velocity about 
the moving center, and the pull on the string is constant. 

(6) A particle is at rest in a smooth horizontal tube. The 
tube is then made to revolve in a horizontal plane with uniform 
angular velocity #. Find the motion of the particle. 

Suggestion. ‘Take the polar coordinates r, ¢, of the particle as 
our coordinates, and let R be the pressure of the particle on 
the tube. 


— =m, 
ar e 
— = mr 
or ? 
oT 

— = mrp 
og 


m[r— ro*] or = 0. 
we (7°h) df = Rrdh. 
dt 
Adding the condition b = at, 


and reducing, ¢—awr=0, 


2morr = Rr. 


12 INTRODUCTION [Ar 6 


Solving, r= Acosh ot + B sinh at, 
r=a and ¢=0 at the start. 
Hence r = a cosh wt = a cosh ¢, 
R= 2 mao’ sinh wt = 2 maw’ sinh ¢. 


If we are interested only in the motions and not in the reac 
tions, problems (a) and (6) can be solved more simply. If in 
each we were to use one less coédrdinate, @ only in (a) and r 
only in (6), rectangular codrdinates X, Y, for the particle could 
be obtained whenever the time was given, and therefore could 
be expressed explicitly in terms of @ or ry and ¢. A careful 
examination of Art. 2 will show that the reasoning is extended 
easily to such a case, and that the work done by the effective 


forces when gq, only is changed is still Feo dg, It is 
1 


to be noted, however, that when the rectangular codrdinates are 
functions of ¢ as well as of q,, ¢,, etc., the energy 7 is no longer 
a homogeneous quadratic in q4,, g,, ete. 


(a') X= nt — asin 6, 
Y=acos@; 
X=n—acos 66, 
Y=—asin 06, 


P= F(X +V = Flv + 0h — 2 an cos 06}. 


= =m [ad — an cos 4), 
oT mere 
39 = mam sin 00. 
m E (a0 — an cos 0) — an sin 06 66 = 0. 
6=0. 


n 
6 =-, as before. 
a 
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(8') T =—[P + 0%"), 
eat = mr 
apd er 
Be = mo’r 
or 


m[7 — wr] dr = 0. 


r =a cosh wt, as before. 


EXAMPLES 


1. A particle rests on a smooth horizontal whirling table ang 
is attached by a string of length a@ to a point fixed in the table 
at a distance 6 from the center. The particle, the point, and 
the center are initially in the same straight line. The table is 
then made to rotate with uniform angular velocity ». Find the 
motion of the particle. 

Suggestion. Take as the single codrdinate @ the angle made 
by the string with the radius of the point. Let X, Y, be the 
rectangular codrdinates of the particle, referred to the line ini- 
tially joining it with the center and to a perpendicular thereto 
through the center as axes. 


Then X= 6 cos ot + a cos(6 + at), 
and Y=06sinot+asin(@ + ot). 
T= 5 [eo + a*(@ + 6)? + 2 abo cos @(@ + 6)], 
o. 2 
and + sind =0; 
and the relative motion on the table is simple pendulum motion, 
the length of the equivalent pendulum being / = . 


2. A particle is attracted toward a fixed point in a horizontal 
whirling table with a force proportional to the distance. It is 
initially at rest at the center. The table is then made to rotate 
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with uniform angular velocity w. Find the path traced on the 
table by the particle. 

Suggestion. ‘Take as coordinates z, y, rectangular codrdinates 
referred to the moving radius of the fixed point as axis of 
abscissas and to the center of the table as origin. Let 4, Y, be 
the rectangular coordinates referred to fixed axes coinciding 
with the initial positions of the moving axes. 


X= x cos ot — y Sin at, Y=-2sin ot + y cos ot. 
T= 5 [2+ we? + y+ wy? — 2 wyt + 2 wry). 
Whence come m[a#— 2 oy — w’x]=— u(a—a), 


m[¥ + 2% — oy] =— py. 


If w? =, the solution is easy and interesting. 
m 


#—2 wy = ao, (1) 
Yt2otr=0. (2) 
Integrating (2), Y¥t2ox=0. 
Substituting in (1), +4 0% = ao’, (3) 


Multiplying (8) by 24, and integrating, 


V+ 4 w’2? = 2 aw’s. 


: jax 
&=20r/)—— 22 
poe 


Whence 2 ot = vers’ ae . 
a 
mao [1— cos 2 wt], 
a 
Veer [2 wt — sin 2 wt]. 
Replacing 2 wt by 0, z= ; [1— cos 6], 


a 
y=—7 [8 —sin 6), 
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and the curve traced on the table is the cycloid generated by a 


circle of radius a rolling backward along the moving axis of Y. 


7. A System of Particles. If instead of a single particle we 
have a system of particles, free, or connected or otherwise con- 
strained, m2, my, mz, are the effective forces on the particle P. 
The effective forces on all the particles are spoken of as the 
effective forces on the system and are mechanically equivalent 
to the set of actual forces on the system. 

T, the kinetic energy of the system, is the sum of the kinetic 
energies of all the particles. 


T= yet +#). 


If SW is the work done by the effective forces in any supposed 
infinitesimal displacement of the system, 


OW = =m [vba + iby + 262). 


If the particles of a moving system are subjected to connec- 
tions or constraints, these connections or constraints may or 
may not vary with the time. In the latter case a set of any n 
independent variables q,, 7,, +++; %,. such that when they and 
the connections and constraints are given, the position of every 
particle of the system is uniquely determined, and such that 
when the positions of all the particles of the system are given, 
%> %** "> Yn follow uniquely, may be taken as codrdinates of 
the system; and m is called the number of degrees of freedom 
of the system. 

In the former case a set of variables q,, 9,,+++, g,) Such that 
when they and the time are given, the position of every particle 
of the system is uniquely determined, and such that when the 
positions of all the particles of the system and the time are 
given, 9,, 9..°**) Y,» follow uniquely, may be taken as the coor- 
dinates of the system; and n is called the number of degrees of 
freedom of the system. 
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The equations expressing the connections and constraints in 
terms of the rectangular codrdinates of the particles and of the 
coordinates 4%, 9. ***) Gs Of the system are often called the 
geometrical equations of the system and may or may not contain 
the time explicitly. In the latter case the geometrical equations 
make it possible to express the coordinates x, y, 2, of every 
point of the system explicitly as functions of the q’s; in the 
former case, as functions of ¢ and the q’s. 

The geometrical equations must not contain explicitly either 
the time derivatives of the rectangular codrdinates of the parti- 
cles or those of the coordinates q,, q,5 +++) qd, of the system 
unless they can be freed from these derivatives by integration. 

Examples of geometrical equations not containing the time 
explicitly are the formulas for transformation of codrdinates in 
Arts. 1 and 8, and the equations for X and Y in Art. 6, (a). 

Geometrical equations containing the time are the equations 
for X and Y in Art. 6, (a’), and in Art. 6, Exs. 1 and 2. 

The work, 6, done by the effective forces when q, is 
changed by 6g, without changing the other q’s is proved to be 


“a or aed 5 
dt eq, eq, za 


by reasoning similar to that used in Art. 2. For the sake of 
variety we take the case where the geometrical equations 
involve the time. 


equal to 


Here e=f [¢, Go Gaels Gu 5 
PROMS Slon LAs BLD 9) ay 4 2 
OS Ope AG aia erates 


and is an explicit function of ¢, @,5 G45 ***s Qn Vy Gort > Ine 





Ot 86008 
- 1 
a, éq, >) 
: d oa ar oe Cin Cin 
and since ——= 





di Og, ~ Bag, + Bg tt ggg, tT gag 


Cuap.I] MOTION OF A SYSTEM OF PARTICLES LT 


02: On Ox . On On 
and = etc laa nig el pam 
dq, 09,06 agi" * ag.aq, > °° + Gq.0q, 
d / ox 0x 2 
zea @) 
8m Ema gh a d 
% °°, 4%, 





fe _a/,e\ ,d/ix\_d ae 
"eq, at Bae Palen al? si) oe Ge 


gov _a[ 2 (#*)_ 0 /# 
“Og, dt|ag,\2/| aq, 5) 


and therefore 6,V= iF Ga = a 84,5 


dt \ eq, oq, 
and if @,q, is the work done by the actual forces when q, is 
changed by 69,, dev eT 
Spe ed Q,- (3) 
0q, 9, 


If the geometrical equations do not contain the time, the same 
result is seen to hold, and in this case it is to be noted that 
since @ is homogeneous of the first degree in the time derivatives 
of the coordinates, that is, in ¢,, ¢,, +++, 9, the kinetic energy 7 
is a homogeneous quadratic in ¢,, 9,,+++) 9. 

Generally every one of the n equations of which equation (8) 
is the type will contain all the n codrdinates 4,, 4, +++, g,, and 
their time derivatives, and can be solved only by aid of the 
others. That is, we shall have n coordinates and the time con- 
nected by m simultaneous differential equations no one of which 
can ordinarily be solved by itself. 

If the forces exerted by the connections and constraints do 
no work, they will not appear in our differential equations. 
Should we care to investigate any of them, we have only to 
suppose the constraints in question removed and the number 
of degrees of freedom correspondingly increased, and then to 
replace the constraints by the forces they exert and to form 
the full set of equations on the new hypothesis. 
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8. A System of Particles. Illustrative Examples. (a) A rough 
plank 16 feet long rests pointing downward on a smooth plane 
inclined at an angle of 380° to the horizon. A dog weighing 
as much as the plank runs down the plank just fast enough 
to keep it from slipping. What is his velocity when he reaches 
its lower end ? 

Here we have two degrees of freedom. Take 2, the distance 
of the upper end of the plank from a fixed horizontal line in 
the plane, and y, the distance of the dog from the upper end 
of the plank, as codrdinates, and let & be the backward force 
exerted by the dog on the plank, and m the weight of the dog. 


ie “ [a +(@+ yy") = = lia" Teh ae, ey). 


‘= =m[22+ 9], 


$ =m[2+y). 


m[2e+y] da = 2 mg sin 80° dx. 
m[a+ i] dy =[R + mg sin 80°] dy. 
2t+Y=49, 


fey ae g 
ae pite Mee 
By hypothesis, x=a constant, 
and therefore PNG? 


2 Hy = 2 WY), 
7 =2Gy+C=2Qy, 
y =V2 gy. 


When y= 16, y = 82, nearly. 
; Rela, 
Since — = 5] 
aa as 2 
_ mg 
Te <a 
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(b) A weight 4m is attached to a string which passes over 
a smooth fixed pulley. The other end of the string is fastened 
to a smooth pulley of weight m, over which passes a second 
string attached to weights m and 2m. 
‘The system starts from rest. Find the 
motion of the weight 4m. 

Two coordinates, z, the distance of 4m 
below the center of the fixed pulley, and 
y, the distance of 2m below the center 
of the movable pulley, will suffice. The 
velocities are ioe aac 
—z for moyable pulley, 
—z+y for 2m, Th 
—az—y for m. 
T=3[4 me + mi*+2m(y—2)?+m(4+y)’) 

=F (8H +37 — 249). 

oT 

= 8z—y], 

ap Lee gl 





oT 

oy 
m[8zé—y]éx=[4 mg — my — 2mg — mg] 82. 
™m [8 y _ r| éy = [2 mg ease mg | by. 





=m[3¥—4]. 


8%—7=0. 
Lode i. 
eB 


The weight 4m will descend with uniform acceleration equal 
to one twenty-third the acceleration of gravity. 

(¢) The dumb-bell problem. Two equal particles, each of mass 
m, connected by a weightless rigid bar of length a, are set 
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moving in any way on a smooth horizontal plane. Find the 
subsequent motion. 

We have three degrees of freedom. Let x, y, be the rectangu- 
lar coordinates of the middle of the bar, and @ the angle made 
by the bar with the axis of X. 

The rectangular coédrdinates of the two particles are 


(2-508 y ~5sin 8) and (2 +5 c0s y+ sin 8); 
their velocities are 
\ | (2 ai sin 06) + (3 = = cos 06) | 
2 2 
and \ | (2 - 5 sin 06) + (u + ; Cos 06) | 
ra") e+ p+ 56 + (asin 0 — a cos) bE+H) 
+P + P+ + (a cos d— asin dd G+ i| 
a . 
aml eee+ Sel 
== 2 m2, 
oT 3 
oy EE 
CL na 
Rin 
2 mxbdr = 0, 
2mydy = 0, 
ma? ~ 
Sp 
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Hence the middle of the bar describes a straight line witb 
uniform velocity, and the bar rotates with uniform angular 
velocity about its moving middle point. 


EXAMPLE 


Two Alpine climbers are roped together. One slips over a 
precipice, dragging the other after him. Find their motion 
while falling. 

Ans. Their center of gravity describes a parabola. The rope 
rotates with uniform angular velocity about their moving center 
of gravity. 


(d) Two equal particles are connected by a string which 
passes through a hole in a smooth horizontal table. The first 
particle is set moving on the table, at right angles with the 
string, with velocity Vag where a is the distance of the particle 
from the hole. The hanging particle is drawn a short distance 
downward and then released. Find approximately the subse- 
quent motion of the suspended particle. 

Let x be the distance of the second particle below its position 
of equilibrium at the time ¢, and @ the angle described about 
the hole in the time ¢ by the first particle. 


T= 5 [e+ a+ (a—2)?6]. 


a i Mz, 
Or 
wee =—m(a— 2) 6, 
0x 
aD =m(a—2)’6. 
06 
m[2¢2+(a—2) 67] da = mgdz, C1) 
m= [(a—2)°6]80 = 0. (2) 
2¢+(a—2)P=9. @) 


(a—2)6=C= aVag, (4) 
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since (2) holds good while the hanging particle is being drawn 
down as well as after it has been released. 


2a4+ ae g = 0, approximately, 


and pd ey 
Qa 


For small oscillations of a simple pendulum of length J, 
d+70=0. 


Therefore the suspended particle will execute small oscilla. 
tions, the length of the equivalent simple pendulum being 2 a. 


EXAMPLE 


A golf ball weighing one ounce and attached to a strong 
string is “teed up” on a large, smooth, horizontal table. The 
string is passed through a hole in the table, 10 feet from the 
ball, and fastened to a hundred-pound weight which rests on a 
prop just below the hole. The ball is then driven horizontally, 
at right angles with the string, with an initial velocity of a 
hundred feet a second, and the prop on which the weight rests 
is knocked away. 

(a2) How high must the table be to prevent the weight from 
falling to the ground? (6) What is the greatest velocity the 
golf ball will acquire? Ans. (a) 8.96 ft. (6) 963.4 ft. per sec. 


9. Rigid Bodies. Two-dimensional Motion. If the particles of 
a system are so connected that they form a rigid body or a 


system of rigid bodies, the reasoning and formulas of Art. 7 
still hold good. 
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(a) Let any rigid body containing a horizontal axis fixed in 
the body and fixed in space move under gravity. Suppose that 
the body cannot slide along the axis. Then the motion is 
obviously rotational, and there is but one degree of freedom. 
Take as the single codrdinate the angle @ made by a plane con- 
taining the axis and the center of gravity of the body with a 
vertical plane through the axis. 

Let 4 be the distance of the center of gravity from the axis, 
and & the radius of gyration of the body about a horizontal 
axis through the center of gravity. Then 


Me a 
vie a (h’ + k*) @. (v. App. A, §§ 5 and 10) 


oT F 
——_- =m hee 6. 
2 ¢ ) 


m (h® + k*) 680 = — mgh sin 086. 
ieee io sin 6 = 0, 
and we have simple pendulum motion, the length of the equiv: 
alent simple pendulum being 
_V+R 
=e eas 

(6) Two equal straight rods are connected by two equal 
strings of length a fastened to the ends of the rods, the whole 
forming a quadrilateral which is then suspended from a hori- 
zontal axis through the middle of the upper rod. The system 
is set moving in a vertical plane. Find the motion. 

Take as codrdinates ¢, the inclination of the upper rod to 
the horizon, and 6, the angle made with the vertical by a line 
joining the point of suspension with the middle of the lower 
rod. From the nature of the connection the rods are always 
parallel. 

Let & be the radius of gyration of each rod about its center 
of gravity. 





U 
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jp eg? a > [a6 ate kd") ae : [2 keg? 4 a6"). 
Gia Neb A) 
ge =2 mi, 
op 


ee = ma’. 
00 
2mk*bdd = 0, 
é=0. 


6+2 sind = 0, 
a 


and the rods revolve with uniform angular velocity while the 
middle point of the lower rod is oscillating as if it were the 
bob of a simple pendulum of length a. 

(ce) If an inclined plane is just rough enough to insure the 
rolling of a homogeneous cylinder, show that a thin hollow 
drum will roll and slip, the rate of slipping at any instant being 
one half the linear velocity. 

Let x be the distance the axis of the cylinder has moved 
down the incline, 6 the angle through which the cylinder has 
rotated, a the radius of the cylinder, and @ the inclination of 
the plane. Call the force of friction /. 


T= 5 (+ #6"). 


ges max. 
ns 
vee — mk?6, 
00 


mise =[mg sin a — F] bz, 
mk?080 = Fadd. 
If there is no slipping, z= a8, 


mi = mg sina — F, 
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mi? = = Fa, 
a 
2 
Hence = = m9 sina — F, 
k? 
te PEE mg Sin & 


2 
For a solid cylinder, k? = a 
F=imgsine; 
R = mg cos @, 
where # is the pressure on the plane; 


Yi a | 
Pee a ees 


3 
where yu is the coefficient of friction. 


For a hollow drum, 


and the drum will slip. 
For the drum, then, 


F =pk = pmg cosa=img sina, 
mz = mg sina —img sina = }mg sina, 
z= %gtsina. 
ma’ = Fa =" mg sin @, 
ab =1gsina, 
af = 1 gt sina. 
S=%—ab = 5 gtsina =, 


where S is the rate of slipping. 
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EXAMPLES 


1. A sphere rotating about a horizontal axis is placed on a 
perfectly rough horizontal plane and rolls along in a straight 
line. Show that after the start friction exerts no force. 


2. A sphere starting from rest moves down a rough inclined 
plane. Find the motion. (a) What must the coefficient of 
friction be to prevent slipping? (4) If there is slipping, what 
is its velocity ? 

Ans. (@) p>#tane. (6) S=gi[sna—fpcosa]. 

3. A wedge of mass M having a smooth face and a perfectly 
rough face, making with each other an angle a, is placed with its 
smooth face on a horizontal table, and a sphere of mass m and 
radius a is placed on the wedge and rolls down. Find the motion. 

Let x be the distance the wedge moves on the table, and y 
the distance the sphere rolls down the plane. 


mla+k ., 2 ea 
9| ge 2 4y cos}. 





t : 
Note that P=5(M+m]e+ 
Ans. (m+ M)x—mycosa=0, Ly—xcosa=4 gt’ sina. 


10. Rigid Bodies. Three-dimensional Motion. (a) A homo- 
geneous sphere is set rolling in any way on a perfectly rough 
horizontal plane. Find the subsequent motion. 

Let 2, y, a, be the coordinates of the center of the sphere 
referred to a set of rectangular axes fixed in space; two of 
which, the axes of X and Y, lie in the given plane. Let O4, 
OB, OC, be rectangular axes fixed in the sphere and passing 
through its center; let OX, OY, OZ, be rectangular axes through 
the center of the sphere parallel to the axes fixed in space; and 
let 6, ¢, be the Euler’s angles (v. App. A, § 8). Take 2, y, 
6, ¢, and y as our codrdinates. The only force we have to 
consider is F, the friction, and we shall let #, and F, be its 
components parallel to the axes OX, OY, respectively. 


rao + YP +POlt el tad), 
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where o,=— Osiny+¢sin 0 cos v, 
Wea 6 cosh + sin @ sin y, 
o, = cos6 +. (v. App. A, § 8) 
Hence P=F[P+ Pt RO+P ty t 2 cos O¢y) J. 
We get me = F., (1) 
my = F,, (2) 
mie = [v + cos 06] = 0, (3) 


mie S (¢ + cos Or] =—aF,sin@siny+aF,sin@ cosy, (4) 
mk? [6 + sin Ody] =— al’, cos p — af, sin 3 (5) 
and as there is no slipping, 
&—aw,=%—a(6 cost ¢sin G sin p)= 0, (6) 
Jtao,=y+a(—Osny+osinOcosp)=0. (7) 
From (4) and (5), 


mk? [sin a (f + cos 04) + sin 4 cosy (6 + sin Oby)] 
=—aF,sn 6, (8) 


mk’ [cos ~ £ (¢ + cos Oy) — sin @ sin y (6 + sin O¢yr)] 
| = af, sin 0. (9) 
Expanding the first members of (8) and (9) and eliminat- 
ing » by the aid of (3), we get 
mk? [cos 6 — sin Pov + sin 6 sin wd + cos 6 sin Lod 
+ sin 0 cos pow] = — aF,, (10) 
mk*[— sin ad — cos WO + sin 4 cos wh + cos 6 cos Wd 
— sin @ sin ver] = aF,. at) 
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But the first members of (10) and (11) are obviously mk’ —* 
Hence, by (6) and (7), 


@,. 


and mk? “° 
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, respectively. 


Substituting in (1) and (2), we get 


whence 


From (7) and (6), 


and from (3), 


Finally, 








mk? .. . 
— ¢ = — mai, 
a 
mh? Se 
- Gi. Y3 
z= 0, 
G0; 
dw 
a 0, 
dt 
dw 
COPE 6 
dt : 
dw 
z—(), 
dt 
B= 0, 
Fo= 0, 
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do, 
dt 


Hence the center of the sphere moves in a straight line with 
uniform velocity, and the sphere rotates with uniform angular 
velocity about an instantaneous axis which does not change its 
direction, and no friction is brought into play after the rolling 


begins. 


(6) The billiard ball. Suppose the horizontal table in (a) is 
imperfectly rough, coefficient of friction p, and suppose the 
ball to slip. 

Take the same coordinates as before, and equations (1), (2), 
(3), (4), (), (8), (9), C10), and (11) still hold good. Let a 
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be the angle the direction of the resultant friction, “= wmg, 
makes with the axis of Y, and let S be the velocity of slip- 
ping, that is, the velocity with which the lowest point of the 
ball moves along the table. Of course the directions of F and S$ 
are opposite. 

Let S, and S, be the components of S parallel to the axes of 
X and Y. We have 

—Scosa=8S,=%—aoa,, 


and —Ssinaw=S,=y% + ao,. 


F, = pmg cos a, 








and F, = pig sin a. 
ds. es as _-« dw 
= —— —=27-—a—. 
7 S sin @& Fie A pg at en 
ds ; da ds dw, 
aa ee ‘7 =Yt+a ar 
From (1), = mg COs a, 
‘dw 
and from (10), Fre =— = Lg COS &. 
2 2 
Hence Ssin aH — cosa = i* Lug COs a, (12) 
and from (2) and (11), 
2 2 
— Scosa@ 46 sing = EE yg sina. (13) 
Multiplying (12) by sin a and (18) by cosa, and subtracting, 
da 
se EN 14 
ss (14) 
Multiplying (12) by cosa and (18) by sina, and adding, 
dS a@+h 
Ebadi (15) 
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2 2 
Integrating (15), S=S,— a — gt. (16) 
From (14), G=,, 


and the direction of slipping does not change. 

If the axes are so chosen that the axis of XY has the direction 
opposite to the direction of slipping, a= 0. Then 

= py, 
Gia): 

These equations are familrac in the theory of projectiles, 
and the path traced on the table is a parabola so long as 
slipping lasts. 

Should y, happen to be zero, the path degenerates into a 
straight line. 

When slipping stops, 


z—ao,=9, and y+ao0,=0, 


and we have the case treated in (a). 


EXAMPLE 


A homogeneous sphere is set rollmg on a perfectly rough 
inclined plane. Find the path traced on the plane. 
Ans. A parabola. 


(ce) The gyroscope. Suppose a rigid body containing a fixed 
point and having two of its moments of inertia about its prin- 
cipal axes through the fixed point equal. Obtain the differ- 
ential equations for its motion under gravity. 

We shall use Euler’s angles with a vertical axis of Z. 


We have o, = 6sin $ — sin 8 cos ¢, 
o, = 6 cos d +> sin @ sin ¢, 

@, = cos 6 + ¢. (v. App. A, § 8) 

T=}[ Ao? + Ao;y + Co?) (v. App. A, § 10) 


=1[4(6 + sin? Of?) + COb cos 6 + 6)?). 
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ae = Cv cos 6 + db), 
ag 
ar ; 
=A sin’ Oy + C cos O(w cos 8+ $), 
oy 
ues — Aé, 
00 
oT 
ie = A sin 6 cos Ov — Csin A (vr cos 8 + dy. 
Our equations are of (+ cos 0+ ¢)= 0, (1) 
. [A sin? Or + C cos 0 (wr cos 8-4 d)]=0, (2) 
A6 — Asin 8 cos OY + C sin 0 (xb cos 6+ $) = mga sin 6. (3) 
From (1), vy cos 6 + d =, (4) 
where @ is the initial velocity about the axis of unequal moment. 
A sin? 0 + Ca cos 0 = L. (5) 
AO — Asin 6 cos OY? + Ca sin OY = mga sin 8, (6) 


or substituting y from (5), 


— Ca cos )(L cos  — Ca) 


(L 
Ab = Asin’? @ 





+ mga sin 6. (7) 


(d) Obtain Euler’s equations for a rigid body containing a 
fixed point. 


Here T =4[Ao? + Bo? + Co;']. (v. App. A, § 10) 
oie Co, (v. App. A, § 8) 
ag 
a = Ao, [6 cos 6 + ap sin @ sin ] 

$ + Bo, [— 6sind + sin A cos] 


= Aw,o, ~ Bo,o,. 
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Whence | ose —(A-—B) 0,0,| do = Nbq, 
¢ 
where NV is the moment of the impressed forces about the C axis. 
dw 
C7 A — B)o,, =—o/Ve 


The remaining two Euler’s equations follow at once from 
this by considerations of symmetry. 


11. In Arts. 2 and 7 it was shown that under slight limi- 
tations the codrdinates of a moving particle or of a moving sys- 
tem could be taken practically at pleasure, and the differential 
equations of motion could be obtained by the application of a 
single formula. It does not follo-v, however, that when it comes 
to solving a concrete problem completely, the choice of coérdi- 
nates is a matter of indifference. Different possible choices 
may lead to differential equations differing greatly in compli 
cation, and as a matter of fact in the illustrative problems of 
the present chapter the codrdinates have been selected with care 
and judgment. That this care, while convenient, is not essential 
may be worth illustrating by a practical example, and we shall 
consider the simple familiar case of a projectile in vacuo. 

Altogether the simplest coordinates are x and y, rectangular 
coordinates referred to a horizontal axis of Y and a vertical 
axis of Y through the point of projection. 


Mtoe 3 
We have fe 5) (+7), 
7 a 
EE es 
ry ce 
mz = 0, 
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Solving, Me 
y=, — gt. 
xr=0,f, 


=F 
y=v,t+ on 


Let us now try a perfectly crazy set of codrdinates, q, and 4, 


where E 
q,=7-+ tan "Ys 


and q, =u —tan-'y. 
Proceeding in our regular way, we have 
=F, + %)s 
2 


hod 


y = tan 
t=3(9,+4,)) 
- PS Eat) Fa A 
y= 5 sec" %1 2 q2 Up aeos 


T= "i+ 24d, + d+ set beh 2g —244,+ +@). 





OL AM 1h—F : 
24, tI 4, +4, +50 BS behG.-a), 
OL ade 4, h— 2 /; -\2 
Se t as > 
ag, 4 se 9 an 9 ng) 


OL ne , QA Oe : 
pa att a — 800 HS 4, - 4), 
m 
4 


el es : 
q ae sepa si cere Clea SAE cere eae AC Hee ‘| 
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tl G+ §,— sect hh een — % tan HOB i i] 


_mg ere 
= Sg POR ctao (2) 


4 


Adding (1) and (2), SC g,) = 9. 
W hence 9, + 4 = 2 Up 

and I+ % = 2 Ute (3) 
Subtracting (2) from (1), 


m 


| sect he died (4, — ¥,) + sec* oii 5 h 2 tan ha _— (4, — %) | 
Ne a9: — 9a. 
=— mg sec mas 
Multiplying by se (4,—4,)» and integrating, 
M0 * 
sect A 5 ‘te (4,—-49,)° =— 8 tan UY 5 24472, 


seo B® G,— 4) = 2 v — 2g tan Le 


it = 
Let 2= arp 
dz eee 
sec ee i 2 9 tan 2, 
ee sec’zdz 


i 
Vv vj — 2g tanz 


1 ae 
t= ly — Vay — 2g tanz] 


i = 
<3] a NE —2gtan h 5 | 





(vy, —gt)* =) } — 2g tan H3s, 
h- 7 t 
an BA = -€, 
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But from (8), tte ane 


9 


a 


Hence G, = Vt + tan7? E ss eI, 


a 


{2 
q, = U,t — tan7} E = a 


Of course this should agree with our first answer, and a 
moment’s consideration shows that it does. 


We have ae =2r=0,, 


2 
tan —o =y=0,t— i , as before. 

12. The parameters UFR eemaor that we have been using to 
fix the position of our moving particle or moving system are 
called generalized codrdinates. Following the analogy of rectan- 
gular codrdinates, the time derivative g, of any generalized 
coordinate g, is called the generalized component of velocity 
corresponding to q, It may be a linear velocity, or an angular 
velocity as in many of our problems, or it may be much more 
complicated than either as in our latest example. 

The kinetic energy 7 expressed in terms of the generalized 
coordinates and the generalized velocities is called the Lagran- 
gian expression for the kinetic energy. 

If we are using rectangular coordinates and dealing with a 


moving particle, he 
T=ZH+H/P+#)3 


oT 


OX 


= Hipitl, 


and is the z component of the momentum of the particle. 


Following this analogy, 2 is called the generalized com- 
qk 
ponent of the momentum of the system, corresponding to the 


coordinate q,. It is frequently represented by p,, and may be a 
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momentum, or a moment of momentum as in many of our prob- 
lems, or it may be much more complicated than either as in 
our latest example. 

Equations of the type 


are practically what are called the Lagrangian equations of 
motion, although strictly speaking the regulation form of the 
Lagrangian equations is a little more compact and will be 
given later, in Chapter IV. 

Q,, defined through the property that Q,8q, is the work done 
by the actual forces when g, is changed by 89,, is called the 
generalized component of force corresponding to q, It may be 
a force, or the moment of a force as in many of our problems, 
or it may be much more complicated than either as in our 
latest example. 


13. Summary of Chapter I. If a moving system has a finite 
number 2 of degrees of freedom (vy. Art. 7) and independent 
generalized codrdinates ¢,, 4,, +++, q,, are chosen, the kinetic 
energy T can be expressed in terms of the codrdinates and 
the generalized velocities q,, Yo. ***y Yay and when so expressed 
will be a quadratic in the velocities, a homogeneous quad- 
ratic if the geometrical equations (v. Art. 7) do not contain the 
time explicitly. 

The work done by the effective forces in a hypothetical infin- 
itesimal displacement of the system due to an infinitesimal 
change 0q, in a single codrdinate g, is 


dor eT), 
dt 04, 09,\ 2" 


Ii this is written equal to @,5g,, the work done by the actual 
forces in the displacement in question, there will result the 
Lagrangian equation 
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The n equations of which this is the type form a set of simul- 
taneous differential equations of the second order, connecting 
the m generalized codrdinates with the time. When the complete 
solution of this set of equations has been obtained, the problem 
of the motion of the system is solved completely. 

It must be kept in mind that in order to obtain the value of 
a single codrdinate or of a set of codrdinates less in number 
than m it is generally necessary to form and to solve the com- 
plete set of m differential equations. 

We shall see, however, in the next chapter, that in certain. 
important classes of problems some of these equations need not 
be formed, and that some of the codrdinates can be safely 
ignored without interfering with our obtaining the values of the 
remaining codrdinates; that, indeed, we may be able to handle 
satisfactorily some problems concerning moving systems haying 
an infinite number of degrees of freedom. 


CHAPTER II 


THE HAMILTONIAN EQUATIONS. ROUTH’S MODIFIED 
LAGRANGIAN EXPRESSION. IGNORATION OF 
COORDINATES 


14. The Hamiltonian Equations. If the geometrical equations 
of the system (v. Art. 7) do not contain the time explicitly 
and the kinetic energy T is therefore a homogeneous quad- 
ratic In g,, Gy +++, Gn» the generalized component velocities, 
Lagrange’s equations can be replaced by a set known as the 
Hamiltonian equations. 

The Lagrangian expression for the kinetic energy we shall 
now represent by 7 


aT; On: 


Let Pia ne P= 5g etc. be the generalized component 
1 2 
momenta. Then p,, p,,+-+ are homogeneous of the first degree 
in 9,,9,.°++. Express ¢,,9,,-+--in terms of p,, p,, +++) 41) or 00s 
noting that they are homogeneous of the first degree in terms 
of p Py +++, and substitute these values for them in 7, which 
will thus become an explicit function of the momenta and the 
coordinates, homogeneous of the second degree in terms of the 
former. This function is called the Hamiltonian expression for 
the kinetic energy, and we shall represent it by Z,. Of course 





T; = T,. () 
By Euler’s Theorem, 
qe 
27;= ait gost t 
4 24, a 24, 
therefore 27,=2T,=p.4,+P.d,+-+* (2) 


eT, 


Let us try to get — and Ep indirectly. 
rim AGS) 
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From (1), 22 = 9054 2%) My, OT Oe 
oq, oq, oq, eq, oq, oq, ; 
ed, | els 00 
or P = — Salt D = oq; + 2 eo 3) 
0g, 09, "189, 8g, ey 
But from (2), 
OT an an 
22a p hy ey... (4) 


aq, dq," P2aq, 
Subtracting (3) from (4), we get 


OTS, OT; 
has me yhapnad! § 5 
24, 29, () 


Again, we have from (1), 
OT, 07; 0g, , OT; Gy 
By 9, op, es P, 








oT OG. n 
or =p - — Reise 6 
op, Pap, ee (6) 
Fr 9 9 OL, oq; OF 
rom (2), eS apd firs Seleay eae (7) 
Subtracting (6) from (7), we get 
Ve 
The Lagrangian equation 
d OT; OT; _ =@ 
dt oq, o% , 
becomes Py + as 0). (9) 
20, 
We have also g aoe (10) 
© ey 


The equations of which (9) and (10) are the type are known 
as the Hamiltonian equations of motion. The so-called canonical 
form of the Hamiltonian equations is somewhat more compact 
and will be given later, in Chapter IV. 
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The 27 equations of which (9) and (10) are the type form 
a system of 27 simultaneous differential equations of the first 
order, connecting the n coérdinates q,, ¢,,+++; q,, and the m com- 
ponent momenta p,, p,,++-+, P~,, With the time, and in order to 
solve for any one coordinate we must generally, as in the case 
of the Lagrangian equations, form and make use of the whole 
set of equations. 

In concrete problems there is usually no advantage in using 
the Hamiltonian forms, but in many theoretical investigations 
they are of importance. It may be noted that in the process of 
forming 7, from 7, q, is expressed in terms of the p’s and q’s, 
and thus equation (10) is anticipated. 

To familiarize the student with the actual working of the 
Hamiltonian forms, we shall apply them to a few problems 
which we have solved already by the Lagrange process. 


15. (a) The equations ot motion in a plane in terms of 
polar coordinates (v. Art. 8, (@)). 


Here i “ [7 + rg). 
ee mr 
Bey es 
oT; A 
—!—p,=mr 
ad Po p 
Whence — Ie, 
pa ay 
femepe 
mr @) 
1 2 
and ——— a Pe O 
OD) a E ot | 
WT aes F 
or mr* 
anae, 
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pdb = Br 8d. 





ats, 
m 
2 
. Pe _ 
Pr mr® : 
d= Pe 
=-*, 
mr 
De Te. 


m[? —rg*] = R, 


our familiar equations. 


41 
(3) 


(4) 


(5) 
(6) 
(7) 
(8) 


(6) Motion of a bead on a horizontal circular wire (v. Art. 5, (@)). 








Here T= x a6", 
a = pia mao, 
06 
pap on, 
mar 
eee 
P29 ma? 


, kp? 
p80 =— ka*6’a80 = — ae 80, 


Donne o: 
pe ma 
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or 











Integrating, ~ 2 ae a =p — 
8 
- maV 
Po Vit +m 
_m Vk 
~ ak Vkt +m 


m kVt 
ees t reise ye ee 
0 ap OS LV* +m]+C= “tog | + Z) 


(¢c) The tractrix problem (v. Art. 6, (a)). 
Here Tes 5 [a? + a’6? — 2a cos Gi}. 


p, = m[& — a cos 06}, 
po = m[a°d — a cos 4). 


; 1 

Whenoe se jaernry: [ap, + cos Op], (1) 
: 1 
= ma? sin? 8 [a cos Op, + po]. (2) 


1 
GP 2a 2 i 
P 9 maz sin? P) [a'p: + Po ae 2 & CO8 Oppo | 


We get p.=sin @, (8) 
: a) 
be [(a’p? + pi) cos 0 + a (1 + cos? 8)p,p9] = 0. (4) 
We have the condition 

ns (5) 
With (5), Pz = m[n — acos 06], 


Pp = ma[ab — n cos 6}. 
Substituting in (4), we get 
Py — mna sin 06 = 0, 
mad = 0, 
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Cuap. IT] 
nv 
Whence Gee es 
a 
Pz, = mn(1— cos 8). 
ed 
bp, = mnsin 06 =" sin 6 = R sin 8. 
a 
2 
mn 
ig mo ’ 
a 
as in Art. 6. 


(d) The problem of the two particles and the table with a 
hole in it (v. Art. 8, (@)). 








Here T= = [2a°+(a—2)**]. 
Pp. = 2M, 
pp =m(a—zx)O 
ny 
M0 Dan. 
Pe 
m(a—x)? 
oes peed 
eA (a—«2)’ 
Pe 
We get SOLE e Bary a 
and De = 0. 
Dye = C= maVag. 
8 
Pot G ol J 
Od, 
Whence ZL era 


if z is small, as in Art. 8, (d). 
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(e) The gyroscope (v. Art. 10, (¢)). 
T, =4[A@ + sin? OY) + Cb cos 6 + $)?)- 


Pp = ut = CC cos 0 + 9), 
op 


py = at = Asin OY + C cos 6 (ar cos 6 + 4), 





ov 
oT; + 
= —!— 46, 
Eee 
Anal ,@ aren zat 
t, =515 [4+ sin? 0 : G 
We get Ps = 9%; @ 
py =9; (2) 
es Asin® 0 [(pé+ py) cos 8 — (14+ cos’ OY py py] = mga sin @. (3) 
Pe C4, 
1 Pia 
pe — Tanta [(C’a® + L”) cos @ — CLa(1+ cos’ 8)] = mga sin 6, 
- Fie (L — Ca cos 6) (L cos 6 — Ca) SGC Tne, (4) 


A sin® @ 
th cos0 +6 =a, 


A sin’? Oy + Cacosd=L, 
as in Art. 10, (¢). 


16. The last two problems have a peculiarity that deserves 
closer examination. Let us consider Art. 15, (e). The kinetic 
energy in the Lagrangian form 7’, and therefore* in the Ham- 
iltonian form 7,, fails to contain the codrdinates @ and wy. 
Moreover, when either of these codrdinates is varied, the 


: Me 

* Since i = ~ <4 (v. Art. 14), it tollows that if a codrdinate is missing in 
Tk Tk 

Tj, it is missing also in T). 
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unpressed forces do no work. Hence two of our Hamiltonian 
equations assume the very simple forms 


Pe =9, Py =9, 
which give immediately 

Ps = Ca, a constant, 
and Py =H, a constant. 


These enable us to eliminate ps and p, from a third Hamilto- 
nian equation (Art. 15, (¢), (3)), which then contains only the 
third codrdinate @ and its corresponding momentum py. 

This same result might have been obtained just as well by 
replacing p, and p, in 7) by their constant values and then 
forming the Hamiltonian equations for @ in the regular way. 
So that if we are interested in @ only, and 7, has once been 
formed and simplified by the substitution of constants for p, 
and p,, the coordinates ¢ and w need play no further part in 
the solution. Should we care to get the values of these ignored 
coordinates, they can be found from the equations p, = Ca, 
py =, by the aid of the value of @ previously determined. 

In Art. 15, (d), since pi, = 0 and py, = ma Vag, the substitution 
of this value for p, in 7, enables us to solve the problem so 
far as x is concerned without paying further attention to 0. 

To generalize, it is easily seen that if the Lagrangian form, 
and therefore the Hamiltonian form, of the kinetic energy fails 
to contain some of the codrdinates of a moving system,* and 
if the impressed forces are such that when any one of these 
coordinates is varied no work is done, the momenta p,, p,, --> 
corresponding to these codrdinates are constant; and that after 
substituting these constants for the momenta in question in the 
Hamiltonian form of the kinetic energy, the coordinates corre- 
sponding to them may be ignored in forming and in solving 
the Hamiltonian equations for the remaining coordinates. 


* Codrdinates that do not appear in the expression for the kinetic energy 
of a moving system are often called cyclic coérdinates. 
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Unfortunately the ignored codrdinates have to be used in 
forming Z;, the Lagrangian form of the kinetic energy, and in 
deducing from it 7, the Hamiltonian form of the energy. 

Not infrequently this preliminary labor may be abridged 
considerably by using a modified form of the Lagrangian 
expression for the kinetic energy of the system, as we shall 
proceed to show. 


17. Routh’s Modified Form of the Lagrangian Expression for 
the Kinetic Energy of a Moving System. In forming the 
Hamiltonian equations of motion (v. Art. 14) we first changed 
the form of 7; by replacing all the generalized velocities 





Gy % ++* by their values in terms of the codrdinates q,, ¢,, 
- and the generalized momenta p,, p,, +++, where Past 
oT; q, 
y=!» ete. 
qo 


Let us now try the experiment of replacing in 7, one only 
of the velocities g, by its value in terms of the corresponding 
momentum p,, the codrdinates q,, g,, +--+, and the remaining 
velocities 9, Gg +++. 

Call 7; thus changed in form, 7,. Of course 

T= Ts, and qd =F (p,, qo» qs9 Poe Cr by 2 *). 

We have OF, _ eT OTs ose 
Tete, seco; 
OL, os oT, pels 








Au ing, —*= oa i,]- 
ransposing mee Prog ~ 89, a i 2 PG) 
aT, aT, 2 
A : (en N —_ Og, 
gain op, oa op, =), ap, ~ bp 5 (ib) dy 
Hence -4=21%,—pa) 
1 op, ‘4 141 


M,, Ta Lo, etal 
is called the Lagrangian expression for the kinetic energy modified 
Jor the codrdinate q,. 
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Our Lagrangian equation 








di'og, 4, 7 
becomes P,- Ee Q: (1) 
oq, 
OM, 
Weh ] 7 = — —4, 
e have also q, a (2) 


It is noteworthy that (1) and (2) differ from the Hamil- 
tonian equations for g, only in that the negative of the 
modified expression M, appears in place of the Hamiltonian 
expression 7, 

Let us go on to the other codrdinates. 

OT, OT; , OT; 2%, 


dpe 


ag, 9, 24, 24,” 








er, eT, 04 a é 
whence Meet le TG eee 
Th Ug A ey eee le 
LEM RUE 
09, 0%, oq, oq, 
ore oT, Wiley eo _, om, 
whence aaa adi ee A eee 
4, Yo 5 Y, Yo 


The Lagrangian equation for g, is therefore 
oe ee Q.» (3) 


and differs from the ordinary form of the Lagrangian equation 
only in that 7, is replaced by the modified expression M,. 

In forming the modified expression it must be noted that 4, 
must be replaced by its value in terms of p,, 9,, Gus Ys 
Yq. ***,) not only in 7, but in the term pg, as well. 

An advantage of the modified form is that when it has once 
been formed we can get by its aid Hamiltonian equations for 
one coordinate and Lagrangian equations for the others. 


48 MODIFIED LAGRANGIAN EXPRESSION [Arr. 18 


The reasoning just given can be extended easily to the case 
where we wish Hamiltonian equations for more than one coérdi- 
nate and Lagrangian equations for the rest. 

The results may be formulated as follows: Let 7,,,y,,-..,», be 
the form assumed by 7; when 4,, 4.+++, G» are replaced by 
their values in terms of p., Py **+> Pro Gr4ar Gr420°* > In Gy 
Gps Gg, Uhen, i 


My, dy 0's yp = Loy dys Pr HE, fe eae — Pps 
we have equations of the type 
ele 
Pp SS a Q de 
é Ou i 
od Gy» Tor 9 Ir 
Up Fa a arn rae cee 
on, 
ifk<rt+1; Pr 
and a O May °° — May ay ste On 
dt oq; 04 
ie JE SS 


18. If we modify the Lagrangian expression for the kinetic 
energy for all the coordinates, 


Mee ca = L SS LeCR, Ie aes SR 
and we get Hamiltonian equations of the form 

iy edi ae, 

fp Oe (1) 
k 0g, & 
: OME wig 

(ee (2) 
‘ OP. 


for all the coordinates, and as we have nowhere assumed in our 
reasoning that 7, is a homogeneous quadratic in the general- 
ized velocities, we can use these equations safely when the 
geometrical equations contain the time explicitly (v. Art. 7). 

If the geometrical equations do not involve the time, in which 
case 7; is a homogeneous quadratic in ¢,, G,) +++, 


2 L; =P, 12.95 ap OO 1p Ge 
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by Euler’s Theorem; and M,,......= 7, —2 7, =— 7,; and (1) 
and (2) assume the familiar forms 
. . aL, 
4 —-=@¢@ 
2 Oly 


It is important to note that the modified Lagrangian expression 


M, _is not usually the kinetic energy of the system, although, 


yy Ign 79 7 
as we shall see later, in some special problems it reduces to the 
kinetic energy. As we have just seen, when the time does 
not enter the geometrical equations, the completely modified 
Lagrangian expression (that is, the Lagrangian expression modi- 


fied for all the coordinates) is the negative of the energy. 


19. As an illustration of the employment of the Hamiltonian 
equations when the geometrical equations contain the time, let 
us take the tractrix problem of Art. 6, (a’). 

Here sie 5 [n? + a°6? — 2 an cos 06] 
and is not homogeneous in 6. 


Do - Tie m[a6 —an cos 9], 





and hes Ps at ge) (1) 
am a 


M,= T, Po Sf ars ae pre ay he 
He = mat sin 8 cos 8+” sin Ope} 
and we have ,— mn’ sin 6 cos 0 — . sin Op, = 0; (2) 


and (1) and (2) are our required Hamiltonian equations. Let 
us solve them. 
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From (1), oe mab — mna cos 6, 
whence pe= mad + mna sin 06. 
Substituting in (2), 
ma?é + mna sin 66 — mn? sin 6 cos 6 — mna sin 66 + mn? sin 6 cos 6 = 0, 
or G= 0, 
which agrees with the result of Art. 6, (a’). 


EXAMPLES 


1. Work Art. 6, (6), by the Hamiltonian method. 
2. Work Exs. 1 and 2, Art. 6, by the Hamiltonian method. 


20. (a) As an example of the employment of the modified 
form, we shall take the tractrix problem of Art. 6 and modify 
for the codrdinate a. 

We have (v. Art. 6) 


(ie 5 [a? + a6? —2acos 6%}. 
Px = m[% — a cos 66}, 


t= P24 a cos 66, 
m 


2 : 
T, = Be + a? sin? 06. 
2 
A= Tl, =p p= |B + a? sin’? 06] _ i —acos 06p, 


m ; ; je e's 
= tle sin’ 06? — Be — 7 @ cos a¢p.|. 








= ma’ sin? 66 — a cos 6p,, 


a) . ° 
a = ma’ sin 6 cos 07 + a sin O6p,. 
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We have for z the Hamiltonian equations 
p, = & sin 0, (1) 
t= P+ a cos 06, (2) 
m 


and for @ the Lagrangian equation 
ma®[sin? 06 + sin 6 cos 06] — a cos 6), = 0. (3) 


Of course (1), (2), and (8) must be solved as simultaneous 
equations, and we can simplify by the aid of the condition 


r= nt. (4) 
Solving, we get 6=0, 
mn 


2 
ra (v. Art. 6 and Art. 15, (¢)) 


(6) As a second example we shall take the problem of the 
two particles and the table with a hole in it (v. Art. 8, (d)) 
and modify for @. 


We have T= = [247+ (a—2x) 6]. 
Po = Ti=m(a —2«)°6, 
whence d= C= (1) 
Pearce 
m= 52+ al - ap 


=Fl22— aig >i 


Our Hamiltonian equations are (1) and 
Po= 0. (2) 
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Our Lagrangian equation is 





lgarcsct 
lee 
aCe) : 
By (2), Pa = C= maVag ; 
whence (3) becomes ae ake oh = Os (4) 
a-—X 


as in Art. 8, (d). 

(c) As a third example we shall take the wedge and sphere 
problem of Ex. 8, Art. 9, and modify for 2. 

We have 


Oe je 





1, = 5 UL+ mye +2] ij 2g cosa]. 

















p.= <= (M+ m) «—my cosa. 
Ld dae MY COS & (1) 
on 
iz pa—my cosa ma+k ., 
Fy 0) (M+ m) is 
M Pa—my cosa ma+kh . p2+mp,y cosa 
z 2(M+m) 2 GP y M+tm ; 


Our Hamiltonian equations are (1) and 


Our Lagrangian equation is 


= mg Sin a, (8) 





E ath m GOs “| . Mp, COS & 





a M+m M+m 
By (2), (ee i6 SV 4 
2 2 A ae} 
whence (8) becomes E _ ie os] ¥ = 9 siD a, (5) 


as in Art. 9, Ex. 3. 
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(4) As a fourth example we shall take the flexible-parallelo- 
gram problem (v. Art. 9, (0)) and modify for ¢. 


We have Hie = [2 kp? ae a6"). 
oT, 9 owe 
= —+ = 2 mk’d. 
Pe a6 p 
ee 
‘ 2 mk? © 
m 2 
eer 3 met | 


Our Hamiltonian equations are (1) and 
Par. (2) 


Our Lagrangian equation is 


mao = — mga sin 6, (3) 
or 6+2sin@=0. (4) 
EXAMPLES 


1. Take the dumb-bell problem of Art. 8, (c), and modify for 8. 





2 
Ans. yy =m [P+ 9 + Pe | 


mas 





2 
M,=m| i? +f oe |: 


mas 
pe = 0. 
mz = 0. 


my = 0. 
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2. Take the dumb-bell problem of Art. 8, (¢), and modify for 
x2 and y. 


m[pitpr oa 
Ans. T,_ »,= ri eres a0 we 
m 24+ pp? 
a Cae aad 
P2= 
Py = 
4 ma = 0. 
3. Take the gyroscope problem, Art. 10, (¢), and modify for 
$ and y. a =f [ses 8 (p —p cos 8)". 
had may 0 C Asin’ @ 


_1| 4g Be Ne Pe 
My = 5) 48 + °C Asin'o |" 


ee i Pepe Mas Be cos 0) cos 0 
Hee C Asin’ @ 


Py = 9, 


Ag — (ps + Pye ss ead 222 ac) Bee Cckaton ee mga sin 0. 
A sin’ 0 

21. We proceed to comment on the problems of the pre- 
ceding section. 

(a) No one of them involves the time explicitly in the 
geometrical equations, and therefore the kinetic energy 7, in 
all of them is a homogeneous quadratic in the generalized 
velocities. 

(6) The momenta, therefore, are homogeneous of the first 
degree in the velocities, and consequently the eliminated veloci- 
ties are homogeneous of the first degree in the corresponding 
momenta and the remaining velocities, and the energy 7,,»,,... 
and the modified function M,,,,,... are homogeneous quadratics 
in the introduced momenta and the velocities not eliminated. 
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(¢) In all the problems the coédrdinates for which we have 
modified the Lagrangian expression for the kinetic energy are 
eyclic. In all of them except the first no work is done when 
any one of the coérdinates in question is varied. Consequently 
one of the Hamiltonian equations for that coordinate is of the 
form p,=0, and the momentum p, = ¢,, where ¢, is a constant. 
Therefore it is easy to express the energy 7,,,,... and the 
modified expression M,,,,... in terms of the remaining codérdi- 
nates, the corresponding velocities, and the constants Ory Cay Pony 
and when so expressed they are quadratics in the velocities but 
not necessarily homogeneous quadratics. When the modified 
function has been so expressed, it may be used in forming the 
Lagrangian equations for the remaining coordinates precisely 
as the Lagrangian expression for the kinetic energy is used, 
and the codrdinates that have been eliminated may be tgnored 
in the rest of the work of solving the problem unless we are 
interested in their values (vy. Art. 16). 

(d) To generalize: If some of the coordinates of a moving sys- 
tem are cyclic, and if the impressed forces are such that when any 
one of these codrdinates is varied no work is done, the momenta 
corresponding to these codrdinates are constant throughout the 
motion. The substitution of these constants for the momenta in 
the Lagrangian expression for the kinetic energy modified for the 
coordinates in question will reduce it to an explicit function of 
the remaining coordinates, the corresponding velocities, and the 
constants substituted, which will be a quadratic in the velocities 
but not necessarily a homogeneous quadratic. 

When the modified function has been so expressed, it may 
be used in forming the Lagrangian equations for the remaining 
coordinates precisely as the Lagrangian expression for the 
kinetic energy is used, and the coordinates that have been 
eliminated may be ignored in the rest of the work of solving 
the problem (vy. Art. 16). 

In the important case where the system starts from rest, the 
constant momenta corresponding to the ignorable codrdinates 
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being zero at the start are zero throughout the motion, and the 
modified expression is identical with the Lagrangian expression 
for the kinetic energy, which therefore is a function of the 
remaining coordinates and the corresponding velocities and is 
a homogeneous quadratic in the velocities (v. Art. 24, (a)). 

(e) The fact that the Lagrangian expression for the kinetic 
energy modified for ignorable codrdinates is expressible in terms 
of the remaining codrdinates and the corresponding velocities 
and is a quadratic in terms of those velocities is often of great 
importance, as we shall see later. 


22. Let us take the gyroscope problem of Art. 10, (c), and 
Art. 20, Ex. 8, and work it from the start, ignoring the cyclic 
coordinates @ and . 

We have 7,=1[4A(@+sin? Oy") + Cop cos6+¢)"], and 
therefore ¢ and y are cyclic codrdinates. Moreover, no work 
is done when ¢ is varied nor when yw is varied, so that ¢ and 
wW are ignorable. 


Pe = 9, and Py =9, so that ewe) and NIE 


We have i= Cap cos6+¢)=e, 
aT, ei : s 
oS sin’ Of + C cos 8 (x cos 0+) = ¢,, 
fi 
" _% _ (4-4 cos 8) cos 8 
pn Reee ? C A sin’ 0 ‘ 
aie Re ee 


A sin? @ 


al gay Cr (Cn IC COMO.) 
T =e By on Gl 2 1 c 
=ere 5| 48 cs C is A sin? @ | 


lft el eC ec0s8) 
My ~ 5| ae 3 
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Forming the Lagrangian equation for @ in the usual way, we 
have 





anes (ce? + ¢F) cos 6 — ee, @! + cos’ 0) 
A sin’ @ 
— ¢, cos 8) (¢, — ¢, cos 0) 
A sin® @ 
which is identical with (7), Art. 10, (¢). 


= mga sin 8, 





or AG + (4 = mga sin 6, 


EXAMPLE 
Work the problem of Art. 20, (6), ignoring the codrdinate 0. 


23. The problem of Art. 20, (@), has an interesting peculiarity. 
Both coordinates are cyclic, and there is no term in the kinetic 
energy that is linear in ¢, the velocity corresponding to the 
ignorable coordinate ¢. Consequently the constant momentum 
Ps 18 a constant multiple of d, which is therefore itself a constant. 
This is true also of p,d and of mk’, the term in the energy 
which involves ¢. T,, and Mf, must then differ from mae” by con- 
stants, and as only the derivatives of J/, are used in forming the 
Lagrangian equation for @, we have merely to disregard the term 
mk? in the energy 7, and use what is left of 7, instead of I/,. 

In cases like this we are able practically to ignore the con- 
tribution to the energy made by the ignored codrdinate as well 
as the coordinate itself, and of course we can conclude that the 
motion we may thus disregard has no effect on the motion we 
are studying, but that the two can go on together without 
interference. 

EXAMPLE 


Examine Exs. 1 and 2, Art. 20, from the point of view of 
the present article. 


24. (a) The wedge and sphere problem of Art. 20, (¢), belongs 
to a very important class. Both x and y are cyclic codrdinates, ana 


x is ignorable. 


Q 


p= t= (+ m)é— my cosa, — p,=0. 
OL 
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The momentum p, is constant, and as it is initially zero 
(since the system starts from rest) it is zero throughout the 
motion, as is p,t. Consequently the kinetic energy 7, and the 
modified expression Jf, are identical, and as they are both 
homogeneous quadratics in y and p, and do not contain y, they 
reduce to the form Ly’, where Z is a constant. 

Therefore our Lagrangian equation for y is Ly =mg sina, 
and the sphere rolls down the wedge with constant acceleration. 
If we care only for the motion of the sphere on the wedge, we 
may then ignore x completely and yet know enough of the 
form of I, to get valuable information as to the required 
motion. Of course we know that the energy of the whole 
system can be expressed in terms of y, and if we are able by 
any means to so express it, we can solve completely for y with- 
out using the ignored coédrdinate w at any stage of the process. 

(6) As a striking example of this complete ignoration of 
coordinates, and of dealing with a moving system having an 
infinite number of degrees of freedom, let us take the motion 
of a homogeneous sphere under gravity in an infinite incom- 
pressible liquid, both sphere and liquid being initially at rest. 

From considerations of symmetry, the position of the sphere 
can be fixed by giving a single codrdinate x, the distance of 
the center of the sphere below a fixed level, and z is clearly a 
cyclic coordinate. 

The positions of the particles of the liquid can be given in 
terms of 2 and a sufficiently large number (practically infinite) 
of coordinates q,, 7,,+++, In a great variety of ways. Assume 
that a set has been chosen such that all the q’s are cyclic.* 
Then, since gravity does no work unless the position of the 
sphere is varied, the q’s are all ignorable. That is, for every one 
of them p, = 0, and the momentum p, = ¢,, and since the system 
starts from rest the initial value of p, is zero, and therefore 
Cea Oe aah , the energy of the system modified for all the g’s, 


CP 


* That this is possible will be shown later, in connection with the treatment 
of Impulsive Forces (v. Chap. III, Art. 36). 


Cuapr. IT] SUMMARY 59 


is then identical with the energy 7,,,,,,... and must be expressi: 

ble in terms of the remaining codrdinate 2 and the corresponding 

velocity 2 (v. Art. 21, (d)); and as @ is cyclic the energy of the 

system will then be of the form Zz’, where Z is a constant. 
Forming the Lagrangian equation for z, we have 


Le = mg, 


and we learn that the sphere will descend with constant 
acceleration. 

Of course this brief solution is incomplete, as it gives no infor- 
mation as to the motion of the particles of the liquid, and since 
we do not know the value of Z, we do not learn the magnitude 
of the acceleration. Still the solution is interesting and valuable. 

The energy of the moving liquid, calculated by the aid of 
hydromechanics, proves to be } m’a*, where m! is one half the 
mass of the liquid displaced by the sphere (Lamb, Hydrome- 
chanics, Art. 91, (3)), and therefore the energy of the system 
is }¢@n+m’)a*, and this agrees with our result. 


25. Summary of Chapter II. The kinetic energy of a mov- 
ing system which has n degrees of freedom can be expressed 
in terms of the n codrdinates q,, q,,+++; J,» and the n general- 
ized momenta p,, P,, +++, Pp» and when so expressed it is a 
homogeneous quadratic in the momenta if the geometrical 
equations do not involve the time (v. Art. 14), and is called 
the Hamiltonian expression for the kinetic energy. If 7, is 
the Hamiltonian expression for the kinetic energy, and 7, the 
Lagrangian expression for the kinetic energy, 





The work done by the effective forces in a hypotheticut 

infinitesimal displacement of the system, due to an infinites). 
3 Meus O, 

mal change 69, in a single coordinate q,, is lit Fee De 
B 

this be written equal to @Q,69,, the work done by the actual 
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forces in the displacement in question, there will result the 
differential equation of the first order, 


: OL, 
Peck 20s = Q,- 


The 2 equations of which this and 
OT, 


> 
op, = es 

given above, are the type are known as the Hamiltonian equa- 
tions of motion for the system. 

If, in the expression 7; — p,4, — 9,9, —+<+— 2.9 Gy Gy °° 9 Gps 
are replaced by their values in terms of p,, Py ++) Po Gav 
=*5 dar Wy > *°*> Gar the result M,,....., is the Lagrangian 
expression for the kinetic energy modified for the codrdinates 
qv ttt) q- For the coordinates for which the expression has 
been modified, that is, for k<r+1, we have Hamiltonian 
equations of the type 


EN gae aera 
;,-— 2M 
: OP, 


For the remaining codrdinates, that is, for k >r, we have 
Lagrangian equations of the type 


doM 0M 
daa Og, elie Cyn Arts £7) 
Whether we work from the Lagrangian expression 7, for the 
kinetic energy, or the Hamiltonian expression 7, or the modi- 
fied Lagrangian expression M,,..., we are in general led to a set 
of simultaneous differential equations whose number depends 
merely upon the number of degrees of freedom in the moving 
system, and such that to solve one we must form and solve all. 
If, however, some of the codrdinates are cyclic (v. Art. 16), 
and the impressed forces are such that when any one of them 
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is varied no work is done, the corresponding momenta are con- 
stant throughout the motion, and these codrdinates are ignor- 
able in the sense that if constants are substituted for the 
corresponding momenta in 7, or M,...., the Hamiltonian equa- 
tions for the remaining coordinates in the former case (or 
the Lagrangian equations in the latter case) can be formed and, 
if capable of solution, can be solved without forming the equa- 
tions corresponding to the ignored codrdinates. 

If the system starts from rest and there are ignorable codér- 
dinates, Jf,,..., the Lagrangian expression modified for the 
ignorable codrdinates, is identical with the kinetic energy of 
the system; and whether the system starts from rest or not, 
M,,,... is a quadratic, but not necessarily a homogeneous quad- 
ratic, in the velocities corresponding to the codrdinates which 
are not ignorable. 


CHAPTER III 
IMPULSIVE FORCES 


26. Virtual Moments. If a hypothetical wfinitesimal displace- 
ment is given to a system, the product of any force by the 
distance its point of application is moved in the direction of 
the force is called the virtual moment of the force, and the sum 
of all the virtual moments is called the virtual moment of the 
set of forces. 

If the forces are finite forces, the virtual moment is the vir- 
tual work; that is, the work which would be done by the forces 
in the assumed displacement. If the forces are impulsive forces, 
the virtual moment is not virtual work but has an interpreta- 
tion as virtual action, which we shall give later when we take 
up what is called the action of a moving system. 


27. For the motion of a particle under impulsive forces we 
have the familiar equations 
M(x, — 2%) = A, 
MY, —Y)= ¥, 
m(é,—4,) = Z. 
Mm(L,— 2), TY =); M (2, — 2), 
are called the effective impulsive forces on the particle and are 


mechanically equivalent to the actual forces. 
If the point is given an infinitesimal displacement, 


m eC, —#,) 6a + Gy.) by + (2,—4,) dz] 


is the virtual moment of the effective forces and of course is 


equal to the virtual moment of the actual forces. 
62 
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If the generalized coérdinates of a moving system acted on 
by impulsive forces are G9 Yo ***, and a displacement caused 
by varying q, by 69, is given to the system, 


8,4 = Em[(é, — 4) ss +(9,—9,) Pay + &,—2 oa 


=Enl G, 4) 2 7 ta yy a, — 4) 5 ma 


where 6, 4 represents the virtual moment of the effective forces. 





Ag im Art. 7, Ee = os 
0, %, 
Therefore ma = me = = a a (#), 
- bea l(ae) — Gn) | 
Hence bale Bah =P, Cl) 


(where P,,.69, is the virtual moment of the impressed impulsive 
forces and P, is called the component of impulse corresponding 
to g,) 1s our Lagrangian equation, and of course we have one 
such equation for every codrdinate 4q,. 

Equation (1) can be written in the equivalent form 


Ce): ck (P),= dS (2) 


28. Illustrative Examples. (a) A lamina of mass m rests on 
a smooth horizontal table and is acted on by an impulsive 
force of magnitude P in the plane of the lamina. Find the 
initial motion. ; 

Let (2, y) be the center of gravity of the lamina, let @ be the 
angle made with the axis of XY by a perpendicular to the line 
of action of the force, and let a be the distance of the force 
from the center of gravity. 
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Then T,=5(@+P +06]. (v. App. A, § 10) 
oT; : 
fey eh a LIES 
Ba er 
P= ey er Ys; 
ai, 24 
= —! = mk’. 
OEE 


Cpe = Cops ae (Poo = 0. 
If the axes are chosen so that 7,=0, y,=0, and 6,=9, 


we have 
me), 


my =F, 
mk’?6 = aP. 
Hence the initial velocities are 


ay ee 0, y = 5 ike 
m mike 


The velocities of a point on the axis of Y at the distance 6 
from the center of gravity are 


V0, 
: : ab| P 
age ib= [144] 5. 
The point in question will have no initial velocity if 6 =— . 
a 


It follows that the lamina begins to rotate about an instan- 
taneous center in the perpendicular from the center of gravity 
2. 
from 





2 
to the line of action of the force at a distance is 


that line and situated on the same side of the line of the force 
as the center of gravity. This point is called the center of 
percussion. 
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(6) A wedge of angle @ and mass MM, smooth below and 
perfectly rough above, rests on a horizontal plane. 

A sphere of radius a@ and mass m is rotating with angular 
velocity Q about a horizontal axis parallel to the edge of the 
wedge and is placed gently on the wedge. Find the initial 
motion (v. Art. 9, Ex. 3). 

Take as codrdinates 2, the distance of the edge of the wedge 
from a fixed axis parallel to it in the horizontal plane; y, the 
distance of the pomt of contact of the sphere down the wedge ; 
and 6, the angle through which the sphere has rotated. 

We have 


Yi 


q 


ze m . 





poe naele 5 — 2ay cosa + k'6’), 


both before and after the sphere is set down. After the sphere 
is set down, y—a0=0. Before it is set down, ¢=y7=0 and 
6=. Since the sphere cannot slip, it exerts an impulsive 
force P up the wedge, and an equal and opposite force P is 
exerted on it by the wedge at the instant the two bodies come 
in contact. 


We have p= ae (M+ m) « — my cos a, 
p=aiem (y — & COS @), 


oT; 28 
= —! = mk’. 
in m6 


Our equations are 








(M+ m) z— my cos a = 0, (1) 

m(y¥ —xcosa)= P, (2) 

mk? (6 —0)=—aP, (3) 

and we have also ab = y. (4) 
From (2), (3), and (4), 

CHE y—teosa= 2. (5) 
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From (1) and (5), 








2 2 
jars m) Z aie ! — mM Cos" ai — i Q cos @, 
p 
2 2 2 
and aren)" sal —m cos? |g = (E+ m) = 0. 
a 


(c) In Art. 8, (6), let the weight 4m be jerked down with a 
velocity v. Find the initial motion. Take the codrdinates x 
and y as in Art. 8, (6), and let P be the magnitude of the jerk. 


T=; BP +37 — 229). 
aT, ae 


7 


ar, iat 
fae me 





Our equations are m(8e—yY)=P, 


and m(sy—«£)=9. 
W hence 23 mz = 3 P,; 
23 my = P. 
But f= Ve e tea 
23 m 
and y= hv 


(ad) Four equal rods freely jointed together in the form of a 
square are at rest on a horizontal table. A blow is struck at 
one corner in the direction of one of the sides. Compare the 
initial velocities of the middle points of the four rods. 

Let m be the mass and 2a the length of a rod. Take a pair 
of rectangular axes in the table. Let (2, y) be the center of the 
figure at any time and @ and ¢ the angles made by two adjacent 
rods with the axis of X 2, y, 0, $, are our generalized codrdinates. 

The rectangular codrdinates of the four middle points are 


obviously (x — @ GOS dh, YSU sin p), qd) 
(v+acos6, y+asin@), Cy 
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(e«+acos¢, y+asin¢), cd) 
We have (w—acos@, y—asin@). (4) 
T,= “ [474474207 +46)+2 (64+ 4%). 
p= 4 md, 
Py = 4 my, 


po =2m(a’ + kh’) 6, 
Py =2me+ k?)d. 
Let the values of 2, y, 0, $, be 0, 0, 9, z, before the blow is 


struck, and let P be the magnitude of the blow. 


Our equations are ami = P, 
4my = 9, 
2m(a+ih)6=0, 

and 2m (a? + k*) b = aP. 
W hence me 
6=0, 

Sg 

4m 

A = ar 


2m(a’ + hk’) 
Let v,, ¥,, U,, ¥,, be the required velocities of the four nuddle 
points. Then 








’ MoO ke ay 
gt tae a +k chin - Sin’ 
Peele De. 

aera ae 8m? 
v,=%—ab= Ege el ei 
8 av+k4m 8 m 
Oy 207 

a ae CUR 


and OE Olt B Oy ORI PA at tsa 
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29. General Theorems. Work done by an Impulse. If a particle 
(x, y, 2) initially at rest is displaced to the position (x + 6z, 
y + oy, 2+ 8z), the displacement can be conceived of as brought 
about in the interval of time é¢ by imposing upon the particle 
a velocity whose components parallel to the axes are u,, 0%, @,, 
where 62 = uot, Oy = vot, and 6z= wot. 

If the particle is initially in motion with a velocity whose 
components are u, v, w, the displacement in question could be 
brought about by imposing upon it an additional velocity whose 





components are obviously u, —u, v, — v, wv, — w. 

Let a moving system be acted on by a set of impulsive 
forces. Let m be the mass of any particle of the system; 
P,, P,, P,, the components of the impulsive force acting on the 
particle; «, v, w, the components of the velocity of the particle 
before, and u,, v,, w, after, the impulsive forces have acted. 

If any infinitesimal displacement is given to the system by 
which the codrdinates of the particle are changed by 62, dy, 
and 6z, we have the virtual moment of the effective forces 
equal to the virtual moment of the actual forces (v. Art. 27): 
that is, 


zm[(u, — u) da + (v, — v) dy + (w, — w) 62] 
=2[P,o¢ + P,dy + P,dz]. 

If the velocity that would have to be imposed upon the parti- 
cle m, were it at rest, to bring about its assumed displacement 
in the time d¢ has the components w,, v,, w,, 2 = uot, dy = v,6t, 
6z=w,6t, and the equation above may be written 

=m [wu 0) 3) 

=2[u,P,+ 0,P, +P, ]. Gb) 


Interesting special cases of (1) are 


zm[(u,—uU)ut+t(, —v)y+(w,—w)w] 
=Z[uP,+vP,+wP,], (2) 
2m[(u, —u)u, +, — 7)», +, — wv), ] 
=> [uP 0 eee) 
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the displacement used in (2) being what the system would 
have had in the time é¢ had the initial motion continued, and 
that in (8) what it has in the actual motion brought about by 
the impulsive forces. If we take half the sum of (2) and (3), 
we get 


Dele +ej+u7]— Ye [aw sLoiyr + w"] 


Ui anes ae » yp ae 
=| See+ P+ tp.) (4) 


2 





or, a systems gain in kinetic energy caused by the action of 
impulsive forces is the sum of the terms obtained by multiplying 
every force by half the sum of the initial and final velocities of 
its point of application, both being resolved in the direction of 
the force. 

This sum is usually called the work done by the impulsive 
forces. 


30. Thomson’s Theorem. If our system starts from rest, 
u=v=w=0, and formulas (1) and (3), Art. 29, reduce re- 
spectively to 

=m [u,u, + ¥,v, + w wJ= ao [(eP +0 wer |, (1) 
and 2m [use of + wy |= 2 (uP, + vP, + f,]. (2) 
But the first member of (1) is identically 


Se {uz + oy +p] + [ul + or + er] 
—[@, Tis aU, my + (w, aN lita 


and subtracting (2) from (1) we get 


m 


9 ate ter] — [alt or + er] 
~[(u,— 4)? + @,— 0)? +(e, — w,)"D 
S(t lente (Dana 0, De (i Ww) Ee | (3) 


If u,, v,, w,, are the components of velocity of the particle m in 
any conceivable motion of the system which could give the 
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points of application of the impulsive forces the same velocities 
that they have in the actual motion, then the second member of 
(3) is zero and we get Thomson's Theorem: 

If a system at rest is set in motion by impulsive forces, its kinetic 
energy is less than in any other motion where the velocities of the 
points of application of the forces in question are the same as in the 
actual motion, by an amount equal to the energy the system would 
have in the motion which, compounded with the actual motion, would 
produce the hypothetical motion. 


31. Bertrand’s Theorem. If Q,, Q,, and Q, are the com- 
ponents of the impulsive force which would have to act on the 
particle m of the system considered in Art. 29 to change its 
component velocities from wu, v, w, to u,, v,, w,, formula (8), 
Art. 29, gives us 


=m [(u, —U) u, + , —v) OPE (w,—w) w,] 
=2[u,Q, + ¥,%, +w,Q,|- (1) 


Subtracting (1) from (1) in Art. 29, we get 


Zm[ (uu, + 2,7, + ww,) — (uy + of + wz)] 
=u, (2-0) ee Ch, = OL) + &, CPs 9) | OQ) 


The first member of (2) is (v. Art. 30) identically 


YF {Led + 0? +002] — [ud + 08 + 2] 
a es es (ee ae mie C0 as w,)*}}. 


If the second member of (2) is zero, as will be the case if the 
Q-forces differ from the P-forces only by the impulsive actions 
and reactions due to the introduction of additional constraints 
which have no virtual moment in the hypothetical motion into 
the original system, we have Bertrand’s Theorem: 

Tf a system in motion is acted on by impulsive forces, the kinetic 
energy of the subsequent motion is greater than it would be if the 
system were subjected to any additional constraints and acted on 
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by the same impulsive forces, by an amount equal to the energy a 
would have in the motion which, compounded with the first motion, 
would give the second.* 


32. By the aid of Thomson’s Theorem many problems in- 
volving impulsive forces can be treated as simple questions in 
maxima and minima. 

(a) If, for example, in Art. 28, (a), instead of giving the 
force P we give the velocity v of the foot of the perpendicular 
from the center of gravity upon the line of the force, so that 
y +0 =», then to find the motion we have only to make the 


OE ca | as . ae 
energy Tr=>[@ +97 +) a minimum. 
aod 


Ges 


a 





ee OPE al : 
de a few+y+ - (v—¥)’]. 


LEO aes) 

6x 

OD ee 8 A : 

—_- = —_— — )— — (0). 

fam 50D] 

Whence reece 

BL ke? F 
bs a+ ih ze 
rem a A 
ce 


and these results agree entirely with the results obtained in 
Art: 28,°(a). 


* Gauss's Principle of Least Constraint: If a constrained system is acted on 
by impulsive forces, Gauss takes as the measure of the ‘' constraint’? what is 
practically the kinetic energy of the motion which, combined with the motion 
that the system would take if all the constraints were removed, would give the 
actual motion. 

It follows easily from Bertrand’s Theorem that this ‘‘constraint’’ is less 
than in any hypothetical motion brought about by introducing additional cons 
straining forces /v. Routh, Elementary Rigid Dynamics, §§ 391-893). 
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(6) In Art. 28, (ec), since =, the energy 
= 


iD = 580+ 389 — 20y). 


To make 7; a minimum, we have 
aT; 
dy 





m[sy—v]=9, 


and the problem is solved. 
(c) In Art. 28, (a), let », = 2+ ab be given. 


Then 7, =F (40, — ab +47 +2@ + PEG 4+$)] 


eo 
oy 
oF; = 2m(a + #)O = 0, 
00 
Ue, eae Bes 5 
—=m[2(¢+k)¢ —4aQ, —4a6)]=9. 
0 
p 4 =0, 
6=0, 
2 av, 
$ B ae ab 
Cae 2 
~ 8e@+h ies age ema 
, pe NS Salles 1 
and Oy EPS ara OF 


33. In using Thomson’s Theorem we may employ any valid 
form in the expression for the energy communicated by the 
impulsive forces. For instance, in the case of any rigid body, 


T= [P+ P+ 2 + do? + Bo} + Co}] 


is permissible and is much simpler than the corresponding form 
in terms of Euler’s codrdinates. 
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Take, for example, the following problem: An elliptic disk 
is at rest. Suddenly one extremity of the major axis and one 
extremity of the minor axis are made to move with velocities 
U and V perpendicular to the plane of the disk. Find the 
motion of the disk. 

Let us take the major axis as the axis of X and the minor 
axis as the axis of Y. 


We have, then, 7 = = [?+e7+te/+ Aol + Bo; + Co;]. 


By the conditions of the problem, since the components of the 
velocity of the point (a, 0, 0) are 0, 0, U, and those of the 
point (0, 6, 0) are 0, 0, V, we have 


z= 0, 
¥yt+ao,= 9, 
z2—aw,= U; 
and & — bo, = 0, 
y= 9, 
2+ bo, = V. 
; BL ee 
Hence r= 3/2 +A +5w—a)], 
or T= wle4 Ola tes oa], 
2 
since a= and B==- 
oem st -)-ju-a|=0, 
s=1U +P). 
1 
We have also 0,=5, 0 Vian), 
1 
c= ho), 
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EXAMPLE 


One extremity of a side of a square lamina is suddenly made 
to move perpendicular to the plane of the lamina with velocity 
U, while the other extremity is made to move in the plane of 
the lamina and perpendicular to the side with velocity /. 


Show that the center will move with velocity 7 perpendicular 


to the plane, and with velocity £ V2 in the plane, toward the 
corner on which the velocity V was impressed. 


34. If our system does not start from rest, it is often easy 
to frame and to solve a problem in which the system is initially 
at rest and is acted on by the same impulsive forces as in the 
actual problem, and where consequently the resulting motion, 
compounded with the actual initial motion, will give the actual 
final motion. 

For example, consider the following problem: A sphere 
rotating about any axis is gently placed on a perfectly rough 
horizontal plane. Find the initial motion. Here, in the actual 
case, the lowest point of the sphere is immediately reduced 
to rest. 

Take rectangular axes of Y and Y, parallel to the plane and 
through the center of the sphere. Let 0,, 0,, ,, be the com- 
ponent angular velocities before, and ,, ,, , after, the sphere 
is placed on the plane. Let 2, ¥, be the velocities of the center 
of the sphere. 

Then, in the actual case, ¢— aw, = 0 and y + aw, = 0 are our 
given conditions. Initially the velocities of the lowest point of 
the sphere are — a, and aQ,. If the sphere were at rest, the 
impulsive force which in the actual case destroys these velocities 
would give to the lowest point the negatives of these velocities ; 
that is, aM, and —aQ,. Let us then solve the following 
auxiliary problem: A sphere is at rest. Suddenly the lowest 
point is made to move with velocities a, and — aQ,,, parallel to 
a pair of horizontal axes. Find the initial motion of the sphere. 
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If w and v are, respectively, the x component and the y 
component of the velocity of the center, and @,, ®,, @,, are the 
angular velocities, 


ie > [ut + v* + kw? + «7 + w;)). 


The velocities of the lowest point are U—aw,, v+ao,, but 
they were given as aQ, and — af... 


Therefore u—ao,=all,, 
v+ao,=— a, 
T= {a* [(@, + 0,)? + (@, + 0,)? + # (@?2 + 02 + @2)]}. 


To make this a minimum, we have 


es =m[a@(o,+ 0,)+ #o,] = 9, 
Cw 


1 


oT , 
Far m [a (@, + 0,,)+ K’o,] = 0, 
‘ ; 
= mk*o, = (0. 
@, - ; 
Hence oO. = g 


— —_——_. (1), 

1 e+ ie x 
2 
a 

0). == = = (0) 

2 a + ke? y 


o, = 0. 


Compounding these with the initial angular velocities in the 


actual problem, 0,, 0,, Q,, we get 
ke? 

Be 
k? 


nt a+ Ke 2,, 

O20): 

These equations, together with ¢—ao,=0 and ¥+ae,=0, 
completely solve the original problem. 
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35. A Problem in Fluid Motion. Let us now consider an inter- 
esting application of the principles of this chapter which was 
made by Lord Kelvin to a problem in fluid motion. 

It is shown in treatises on hydromechanics that if an incom- 
pressible, frictionless, homogeneous liquid, either infinite in 
extent or bounded by any finite closed surfaces fixed or moving, 
and with any rigid or flexible bodies immersed in it, is moving 
under the action of conservative forces (v. Chap. IV) and has 
ever been at rest, the motion will be what is called irrotational. 
That is, if x, y, z, are the rectangular coordinates of any fixed 
point in the space occupied by the liquid, there will be a func- 
tion ?(2, y, z) such that if u, v, w, are the components of the 
velocity of the liquid at the point 2, y, 2, 

é 0 0 
=~, = s wa. 
The function ¢ is called the velocity-potential function. 

Since throughout the motion the liquid is always supposed 
to be an incompressible continuum, u, v, w, must satisfy the equa- 
tion of continuity for an incompressible liquid, 





du ov, ow _g, 
Ox Oy Oz 
and therefore ¢ satisfies Laplace’s equation, 
ap Fh od 
Ox” + by? a ee: 


and will be uniquely determined except for an arbitrary constant 

‘ 0 : , 

term if the value of —, the velocity normal to the surface, is 
on 


given at every point of the boundary of the liquid, however 
irregular that boundary. Therefore the actual motion at every 
point of the liquid at any instant is uniquely determined, if the 
motion is irrotational, when the normal velocities at all points 
of the boundary are given. 

We wish now to prove that the kinetic energy of the actual 
motion is less than that of any other motion, not necessarily 
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aTotational, consistent with the equation of continuity and with 
the actual normal velocities at the boundary. 

If u, v, w, are the velocities at (a, y, 2) in the actual motion, 
let uta, v+ 8, w+y, be the velocities in the hypothetical 
motion, and let v, be the actual normal velocity at any point 


$ 


6 
of the boundary. Then we have == lw +mv+nw =v,, where 
on 


l, m, n, are the direction cosines of the normal, and wu, v, w, are 
the components of the velocity at the point in question. 
In the hypothetical motion, 
l(uta)+mv+B8)+nwt+y)=%, 
at the same point; therefore 
la+mB+ny=0 (1) 
at every point of the boundary. 
As the hypothetical velocities as well as the actual velocities 
must obey the law of continuity, 
é(ut+ea) eo AEs B) ea 6 Ss Woe 0, 
ox 0. 02 
er 08 | Oy 
CORR Tol) 2) 2 
and therefore ae + ope ae (2) 


at every point in the bounded space. 
If 7 is the energy of the actual motion, 


rab [fire te' dx dy daz, 


where p is the density of the liquid, and where the volume 
integral is taken throughout the space filled by the liquid. 
If 7’ is the energy of the hypothetical motion, 


r= 8 ([[lutart@o+ By +(w+)'] dedy de 


= 248 [[ [late +9'] dedy de 
+e fff laut Bet re] dx dy dz. 
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By the aid of Green’s Theorem we can prove that 


il [au + Bu + yw] dxdydz= 0. 


We have We dx dy dz = [ Veos ads 
xv 


(v. Peirce, Newtonian Potential Function, p. 92 (143)), where 
the volume integral is taken throughout any bounded space 
and the surface integral over the boundary of the space, cos a 
being the 2 direction cosine of the normal to the boundary. 


Now Ge ee a 


and ff fe dx dy dz = | lagas ain — de dy de. 
Ay 


In like manner, 


f[ffeedeayae= [ magas— [fot dedyae 
oe [ff redeay de = nvaas— [fH dvdy ae 


Hence iif [aw + Bu + yw] da dy dz 


= [et me-+mngas— ff fo] E42 +2] dedyae 


But the surface integral vanishes by (1), and the volume 
integral vanishes by (2). Therefore the energy Z’ is greater 
than the actual energy 7. 

It follows that the irrotational motion of any frictionless 
ncompressible homogeneous liquid under the action of con- 
servative forces is at every instant identical with the motion 
which would have been suddenly generated from rest by a set of 
‘mpulsive forces applied at points in the boundary of the liquid 
and such that they would suddenly give all the points of the 
boundary the normal velocities that these points actually have 
at the instant in question (v. Thomson’s Theorem in Art. 30). 
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36. If now we have a liquid contained in a material vessel 
and containing immersed bodies, a set of generalized codrdinates 
I> I °° > Guy CaN be chosen, equal to the number of degrees of 
freedom of the material system formed by the vessel and the 
immersed bodies, and the normal velocity of every point of the 
surface of the vessel and of the surfaces of the immersed bodies 
can be expressed in terms of the codrdinates q,, 9, +++, g,, and 
the corresponding generalized velocities ¢, g++ +s Gy 

We can now choose other independent codrdinates gi, gj, +++, 
practically infinite in number, which, together with our codrdi- 
nates 9, 9.°**, G,» Will give the positions of all the particles 
ot the liquid. 

Suppose the system (vessel, immersed bodies, and liquid) at 
rest. Apply any set of impulsive forces, not greater in num- 
ber than n, at points in the surface of vessel and of immersed 
bodies and consider the equations of motion. For any of our 
coordinates qg, we have the equation 


P= 9, 


where p,; is the generalized momentum corresponding to gj, 
since (as in varying gq no one of the coordinates ¢,, ¢,, +++, dy 
is changed, and therefore no one of the impulsive forces has its 
point of application moved) the virtual moment of the com- 
ponent impulse corresponding to q is zero. 

As the actual motion at every instant could have been gen- 
erated suddenly from rest by such impulsive forces as we have 
just considered, the momentum p, is zero throughout the actual 
motion; and the impressed forces being by hypothesis consery- 
ative, and the liquid always forming a continuum, no work is 
done when gq, is varied. Consequently, in the Hamiltonian 





equation p, + ae =, p,=90 and Q,= 0, and therefore a = 0. 
k 


Hence every coérdinate 9g, is cyclic, and it is also completely 
ignorable. The energy modified for the coordinates gq; is identi- 
cal with the energy, which, being free from the coordinates q, 
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and the momenta p,, is expressible in terms of the coordinates 
dy Yq °°) Gy and the corresponding velocities g,, G5 +++ du 
and is a homogeneous quadratic in terms of these velocities 
CVs Artec (a).)s 


37. Summary of Chapter III. In dealing with problems in 
which a moving system is supposed to be acted on by impul- 
sive forces, we care only for the state of motion brought about 
by the forces in question, since on the usual assumption that 
there is no change in configuration during the action of the 
impulsive forces, we are not concerned with the values of the 
coordinates but merely with the values of their time derivatives. 

The virtual moment (v. Art. 26) of the effective impulsive 
forces in a hypothetical infinitesimal displacement of the system, 
due to an infinitesimal change 64, in a single codrdinate q,, is 


i Al Ea ) 5 
oa. | lag | pode OF [Peds —CPedol 04. 
lalla} & (ee. -@oae 
If either of these is written equal to P,8q,, the virtual moment 
of the actual impulsive forces in the displacement in question, 





we have one of the equivalent equations 


lal Lea) =? 
Ch}, Led do 
Chk= (Pr) =P. 

The n equations of which either of these is the type are » 
simultaneous linear equations in the x final velocities (4,),, 
(4,) ***, and as the configuration and the initial state of 
motion are supposed to be given, no integration is required, 
and the problem becomes one in elementary algebra. 

A skillful use of Thomson’s or of Bertrand’s Theorem re- 
duces many problems in motion under impulsive forces to 
simple problems in maxima and mimima. 


CHAPTER IV 
CONSERVATIVE FORCES 


38. If X, Y, 4, are the components of the forces acting 
on a moving particle (coédrdinates 2, y, 2), the work, W, done 
by the forces while the particle moves from a given position 
Py, (@y Yor %)» to a second position P,, (2,, ¥,, 2,), is equal to 


Py 
P [Xdx + Ydy + Zdz]; 
Po 


and since every one of the quantities X, Y, and Z is in the 
general case a function of the three variables z, y, z, we need 
to know the path followed by the moving particle, in order to 
find W.. Let fG@, y, 2=0, 6@ y, 2) = 9, be the equations: of 
the path. We can eliminate z between these equations and then 
express y explicitly in terms of x; we can then eliminate y 
between the same two equations and express z in terms of 2; 
and we can substitute these values for y and z in X, which will 


wy 
then be a function of the single variable 2, and i Xdx can be 
TM, 
found by a simple quadrature. By proceeding in the same way 


with Y and Z we can find {_ YVdy and fZde, and the sum 


Ny 

of these three integrals will be the work required. 
It may happen, however, that XYdx + Ydy + Zdz is what is 
ealled an exact differential, that is, that there is a function 


U=$¢(@, y, 2) such that 
6x oy 02 
81 
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Since the complete differential of this function is 


0U oU 
Ox oy 


or Xdx+ Ydy + Zdz, we have 


(xa + Ydy + Zdz] =o, y, D=U, 


oU 
“= dz, 
Mpa Z 


and apace + Ydy + Zdz] 

an = $y I 4%) - $ (Xs Yope Na U5 
and in obtaining this result we have made no use of the path 
followed by the moving particle. 

When the forces are such that the function U=4(%, y, 2) 
exists, they are said to be conservative, and U is called the 
force function. 

We can infer, then, that the work done by conservative 
forces on a particle moving by any path from a given initial 
position to a given final position is independent of the path and 
is equal to the value of the force function in the final position 
minus its value in the initial position. 

If instead of a moving particle we have a system of particles, 
the reasoning given above applies. 

Let (2%, Y» %) be any particle of the system, and Y,, Y,, Z,, 
be the forces applied at the particle. Then the whole work, W, 
done on the system as it moves from one configuration to 
another is equal to 


(Pe), 
s iL [X,du, + Y,,dy, + Z,dz,]. 


Bo 


If there is a function V=¢$(a, @, +++) tyres Yy Yo ees 
Ye ez zy Zo ene asi am) oe -) such that 
Cus CU. OU 
ox, > Een? OY, Tek? 02, Tine es 


then 2 [X,dx,+ Y,dy,+Z,dz,] is an exact differential and 
U =f (By By 0 ey Yur Yor 0 y By Z,,+++) 1s its indefinite integral 
and is the force function; the forces are a conservative set; and 
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the work done by the forces as the system moves from a given 
initial to a given final configuration is equal to the value of 
the force function in the final configuration minus its value in 
the initial configuration, no matter by what paths the particles 
may have moved from their initial to their final positions. 

It is well known and can be shown without difficulty that 
such forces as gravity, the attraction of gravitation, any mutual 
attraction or repulsion between particles of a system which for 
every pair of particles acts in the line joining the particles and 
is a function of their distance apart, are conservative ; while such 
forces as friction, or the resistance of the air or of a liquid to 
the motion of a set of particles, are not conservative. 

The negative of the force function of a system moving under 
conservative forces is called the potential energy of the system, 
and we shall represent it by V. 

If we are dealing with motion under conservative forces and 
are using generalized codrdinates, and the geometrical equations 
do not contain the time, we can replace the rectangular coordi- 
nates of the separate particles of the system in the force function 
or in the potential energy by their values in terms of the gener- 
alized coordinates 9,, ¢,,++-, g,, and we can thus get U, and con- 
sequently V, expressed in terms of the generalized codrdinates. 
If U is thus expressed, Q,69, (the work done by the impressed 
forces when the system is displaced by changing 4g, by 6q,) is 6,,U 


Tk 
‘ : oU 0 

and therefore is approximately i, be or = dq, and hence 

k k 


OP 


39. The Lagrangian and the Hamiltonian Functions. If the 
forces are conservative, our Lagrangian equation 


C (1) 








wri — = ——, 2 
may be written diog, 08. 0 OG, (2) 
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where V is the potential energy expressed in terms of the coordi- 
nates gy G++ *s and not containing the velocities ¢,, ¢,,-+- >. 

If dS ad (3) 


L is an explicit function of the codrdinates and the velocities 
and is called the Layrangian function. 





Obviously Ub an ole aaa a z on om g Ve 
OG, 9% OG, OGyn 1003 


Hence our Lagrangian equation (1) can be written very neatly 


a foil OL 
— | —_- j= —_.- 4 
dt bal O”% D 


If our forces are conservative, and we are using the Hamil- 
tonian equations and express the kinetic energy 7, in terms of 
the codrdinates and the corresponding momenta, and if we let 


H=T,+YV, (5) 
H is called the Hamiltonian function and is a function of 4, 4,, 


=*) Jy and Dra Danco Dee 
Our Hamiltonian equations 
CL, 
Pet a = VO, 


k 


one 
Cp e= OD, ’ 
can now be written 
dp, Ke CHEE 
ih oq, ©) 
dq, 0H 
whe Cp, <> 


and these are known as the Hamiltonian canonical equations. 

If our forces are conservative, and we are using instead of 
the kinetic energy the Lagrangian expression for the energy 
modified for some of the codrdinates q,, g,,+++) g (Vv. Art. 17), 


d if ] 
and ut we let Cy aes 5 (8) 
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we have for any coordinate g,,, the Lagrangian equation 


(22). #2 oS 
dt 0”, +3: OG, 4% ; 


and for any codrdinate g,_, the pair of Hamiltonian equations 








Ap, , Oo 

es (10) 
AG, p = o® 

dt a OP,-- oe 


40. 'The Lagrangian function Z is the difference between the 
kinetic energy 7; (expressed in terms of the coérdinates 4, ¢,, 

-+, g,, and the velocities q,, ¢,, +++; G,, and homogeneous of the 
second degree in terms of the velocities) and the potential 
energy V (expressed in terms of the codrdinates 4,, 9,,--+, ¢,)- 

The Hamiltonian function H is the sum of the kinetic 
energy 7, (expressed in terms of the codrdinates q,, 9, +++ dus 
and the corresponding momenta p,, p,,--+, p,, and homogen- 
eous of the second degree in terms of the momenta) and 
the potential energy V (expressed in terms of the coordinates 
> Yor **% Gn) 

The sum of the kinetic energy and the potential energy, 
however expressed, is sometimes called the total energy of the 
system, and we shall represent it by #, so that 


Biz! VY. (1) 
The function ® of the preceding section is the difference 
between the kinetic energy (expressed in terms of the codrdi- 


nates 9,5 9%. °**> Gar the momenta p,, Py +++, Py» and the 
velocities 9,415 G42. °**» %,)» minus the terms 


Pits t+ Pale +++ + Prd 
(similarly expressed), and the potential energy (expressed in 


terms of the codrdinates q,, 9,, +++» q,)- We shall call it the 
modified Lagrangian function. 
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It is to be observed that all the terms of ® except those 
contributed by the potential energy are homogeneous of the 
second degree in the momenta introduced and the velocities not 
eliminated by the modification. 


41. In dealing with the motion of a system under conserva- 
tive forces, we may form the differential equations of motion 
in any one of three ways, and the equations in question are 
practically given by giving a single function — Z, the Lagrangian 
function, or H, the Hamiltonian function, or ®, the modified 
Lagrangian function.* 

Every one of these functions consists of two very different 
parts: one, the potential energy V, which depends merely on the 
forces, which in turn depend solely on the configuration of the 
system; the other, the kinetic energy Z or the modified 
Lagrangian expression M, ..., either of which involves the veloci- 
ties or the momenta of the system as well as its configuration. 

If we are using as many independent codrdinates as there are 
degrees of freedom, a mere inspection of the given function 
will enable us to distinguish between the two functions of 
which it is formed, the potential energy or its negative being 
composed of all the terms not involving velocities or momenta. 

If, however, we are ignoring some of the coordinates (v. Arts. 
16, 21, and 24) and are using H (the Hamiltonian function) or ® 
(the Lagrangian function modified for the ignored codrdinates), 
the portion contributed to H by 7, or to ® by M (the modified 
expression for the kinetic energy) is no longer necessarily a 
homogeneous quadratic in the velocities and momenta (v. 
Art. 21) and may contain terms involving merely the codrdi- 
nates and therefore indistinguishable from terms belonging to 
the potential energy; and consequently the part of the motion 
not ignored would be identical with that which would be pro- 
duced by a set of forces quite different from the actual forces. 


* Indeed, for equations of the Lagrangian type, any constant multiple of L 
or @ will serve as well as L or ®, 
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We may note that the last paragraph does not apply if the 
system starts from rest, so that the ignored momenta are zerc 
throughout the motion (v. Art. 24, (a)). 

Let us consider the problem of Art. 8, (d), where the potential 
energy V is easily seen to be — mgz. 

If we are using the Lagrangian method, as in Art. 8, (a), 
we have 


i ae ‘ 
L=G[2e+(a—a)@ + 2gz). ER 
If we use the Hamiltonian method, as in Art. 15, (d), we 
have | 
Px pe 2 
H= 2. 2 i 2 
2 ale + (a—«x)’ me a] 2 


If we use the Lagrangian method modified for 6, as in 
Art. 20, (6), we have 


We De 
® = —| 2¢?-—_“* __+ 2 gz}. 3 
3 ¢ ay Cran yas oe] Sy 
A mere inspection of any one of these three functions enables 
us to pick out the potential energy as V =— mgz. 
If, however, we aré ignoring 0, as in Art. 22, Example, we 
have 2 
ee 2 oe 4 
m (a Se J | SS. 


and here, so far as the function ® shows, the potential energy 
2 


may be — mgz, as, in fact, it is, or it may be Mga. 


nti ay i * 
Indeed, the hanging particle moves precisely as if its mass were 
2m and it were acted on by a force having the force function 


Ce 
—————— MNGXLs 
Im(a—zyi* cag 
2 
m(a—x)* 
If, as in many important problems (v. Chap. V), we are unable 
to discern and measure the impressed forces directly and are 
attempting to deduce them from observations on the behavior of 


that is, a force vertically downward equal to mg — 
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a complicated system, which for aught we know may contain 
undetected moving masses, the fact that we cannot discriminate 
with certainty between the terms contributed to the modified 
Lagrangian function by the kinetic energy of the system and 
those contributed by the potential energy may lead to entirely 
different and equally plausible explanations of the observed 
phenomena (v. Art. 51). 


42. Conservation of Energy. If we are dealing with a system 
moving under conservative forces, the coordinates 9,, 9,, +++) Yu» 
are functions of ¢, the time, as are also the velocities 4,, 9,, -- +5 4,- 
Therefore V, the potential energy, and 7, the kinetic energy, are 
functions of the time, as is their sum, the Hamiltonian function H. 


Let us find es 
dt 


As H depends explicitly on the codrdinates q,, q,, +++, G,. and 
the momenta p,, P.5 +++) Pus 
dH OH dq, _ 0H dp, 
pipes : if 
dt lhe dt i, dt > 
But by our Hamiltonian equations (v. Art. 389, (6) and (7)), 
dp, __ eH dy, _ eH 


— » and 5 
dt On, di op, 








and (1) becomes 





aH Olsl Gisl OE IEE 

sere a = 2 

dt a Op, OP, Al SS 
Therefore T+V=H=h, (3) 


where / is a constant. 

Hence in any system moving under conservative forces the 
sum of the kinetic energy and the potential energy ts constant 
during the motwon. 

This is called the Principle of the Conservation of Energy. 

Since by (3) any loss in potential energy during the motion 
is Just balanced by an increase in the kinetic energy, and the 
loss in potential energy is equal to the work done by the actual 
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forces during the motion, our principle is a narrower statement 
of the familiar principle: Jf a system is moving under any forces, 
conservative or not, the gain in kinetic energy is always equal to 
the work done by the actual forces. 


43. Hamilton’s Principle. Let a system move under conserv- 
ative forces from its configuration at the time t, to its con- 
figuration at the time ¢. We have 

eae (1) 
dt Od, FO, 
where ZL, the Lagrangian function, is equal to 7— V. 

Suppose that the system had been made to move from the 
first to the second configuration so that the particles traced 
slightly different paths with slightly different velocities, but so 
that at any time ¢ every codrdinate q, differed from its value 
in the actual motion by an infinitesimal amount, and so that 
every velocity g, differed from its value in the actual motion 
by an infinitesimal amount, or (using the notation of the 
calculus of variations)* so that 6g, and 64%, were infinitesimal ; 
and suppose it had-reached the second configuration at the 
time ¢,. Then, if at the time ¢ the difference between the value 
of Z in the hypothetical motion and its value in the actual 
motion is 6, 


OL OL 
oL = — 69, + — oq | 2 
> Pera 2) 
Now, at the time ¢, 
ip eee 
Tk dt Tk 
anil ey, eC 


aj, Gq, dt 
d foL d oL 
= oq, é 
agi i) dt oj, 


d /oL oL 
= <= a é 5) by 1 © 
Ae a) 0, G by CA) 


* For a brief introduction to the calculus of variations, see Appendix B. 
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d OL 
= —_— —— 8 9 
Therefore, by (2), 6Z a > EF | 


4 t t 
and if Ves if rat =| 37-5, | 
ty ty k ty 


Since the terminal configurations are the same in the actual 
motion and in the hypothetical motion and the time of transit 
is the same, 6g,=0 when ¢=#, and when ¢=1¢, and 


3 { 1ae=df [r-v]ae=0. (3) 
t 7 


But (8) is the necessary condition that Z should be either a 
maximum or a minimum and is sometimes stated as follows: 
When a system ts moving under conservative forces, the time inte- 
gral of the difference between the kinetic energy and the potential 
energy of the system is “ stationary” in the actual motion. This 
is known as Hamilton’s principle.* 


44, The Principle of Least Action. If the limitation in the pre- 
ceding section that “in the actual and the hypothetical motions 
the time of transit from the first to the second configuration is 


* Hamilton’s principle plays so important a part in mechanics and physics 
that it seems worth while to obtain a formula for it which is not restricted to 
conservative systems. 

We shall use rectangular codrdinates, and we shall make the hypothesis as to 
the actual motion and the hypothetical motion which has been employed above. 

For every particle of the system we have the familiar equations 


(ite = 2G, Wi = Ve bs 3 


Since T= » . (#2 + 92 + 2%), 





dT = Em [hse + dy + 262] = zm| 2 de + i Soy + 202} 
But mi a 62 = mM £ (£6x) — mada = f (madx) — Xda. 
dt dt at 
d 
dt 
and i “(974 B[Xoe + Yby + Zde]} dt = B [mbox + poy + Zaz} = 0. 
0 


Hence 67 + 2[X6e + Yéoy + Zéz] = — Xm [sda + Ydsy + 282], 


If the system is conservative, 2[ X62 + Ydy + Zz] = 5U = — SV, and we 
get the formula (8) in the text. 
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the same” be removed and the variations be taken not at the 
same time but at arbitrarily corresponding times, ¢t will no longer 
be the independent variable but will be regarded as depending 
upon some independent parameter 7. We now have (v. Art. 43) 


d d 
B3 = . App. B, § 6, (1 
oo dt og, — dk dt ot, (v. App. B § 6 (1)) 
aLy. hd,  .obds, 
ag, 1 ag, at 2 296, dt 
d /eL aL aL d 
Stef OG Hy |e eae Oe een OF. 
le a) oq, 1* 104, dt” 
d <a 1eL d aL 
ee RE ee ey ee ce 
ae a als | di > |é 7 
dice ROL 
ue Oe ea NG nN. 
ala, | d 
OL OT, 
since a= as 
Ode =6P"FU 
d Mae 
pes Eira 
Hence eee a, t olen | 


: Bea pec ie 
and ip jez tore eee i errs Sain 
OL t 
=| D552 =0. eb) 


If now we impose the condition that during the hypothetical] 
motion, as during the actual motion, the equation of the con 
servation of energy holds good, that is, that 


TEV Sh, 
then 1) 0) — 9, 
8L = 87 —8V =287 =8(2T), 
and (1) becomes 
[Panter Z| a=0 
| dt ld 
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if fe art +27 w]ar=0, 


{OE SQT)+2 79S |dr= 0, 
JiAe5)* Py 
ef one dr = 0, 


and, finally, 3 f° Data); (2) 


a 
The equation A =/ 27 dt defines the action, A, of the forces, 


ty 


and the fact stated in (2) is usually called the principle of least 
action. As a matter of fact, (2) shows merely that the action is 
“ stationary.” 


45. In establishing Hamilton’s principle we have supposed 
the course followed by the system to be varied, subject merely 
to the limitation that the time of transit from imitial to final 
configuration should be unaltered; consequently, as the total 
energy is not conserved, the varied course is not a natural 
course. That is, to compel the system to follow such a course 
we should have to introduce additional forces that would 
do work. 

In establishing the principle of least action, however, we 
have supposed the course followed by the system to be varied, 
subject to the limitation that the total energy should be unal- 
tered; consequently the varied course is a natural course. That 
is, to compel the system to follow such a course we need 
introduce merely suitable constraints that would do no work. 

We have deduced both principles from the equations of 
motion. Conversely, the equations of motion can be deduced 
from either of them. Each of them is therefore a necessary 
and sufficient condition for the equations of motion. 
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46. The action, A, of a conservative set of forces acting on a 
moving system has been defined as the time integral of twice 


the kinetic energy. ; 
De i 2 Tilt. (1) 
t 


0 


% ty o\ 2 t 
But A =[ Im(e?+y?+#)dt =/ Sm (=) at -[ Smv7dt. 
t ty t 


0 


Therefore A = [ Smnds (2) 


so that the action might just as well have been defined as the 
sum of terms any one of which is the line integral of the momen- 
tum of a particle taken along the actual path of the particle. 

There is another interesting expression for the action, which 
does not involve the time even through a velocity. 


Smee T+V=A, 27=2(h—-YV). 
As T=>>@+H+®), 
ime +y+#)=2(A—-V), 
1s\? 
Sl | =2h— 7), 
n( 5) (h ) 
=mds? = 2(h— V) dt’, 
| =mds* 
and dt= op ay 
=mds* 
= | 2(h-V —_—_—__. 
Behe esON sry) 


4 ji J2(h— Vy Smads. (3) 





47. We have stated without proof that the differential equa- 
tions of motion for any system under conservative forces could 
be deduced from Hamilton’s principle or from the principle of 
stationary action. Instead of giving the proof in general, we will 
give it in a concrete case, that of a projectile under gravity. 
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We shall use fixed rectangular codrdinates, Y horizontal and 
X vertically downward, taking the origin at the starting point 
of the projectile. In this case, obviously, 


T= = (@+%), and V=— mga. 
(a) By Hamilton’s principle, we have 


Be sh key : 
ef sl? +e + 292] d= 0, @) 
0 a 
or f Set 9 +292) d= 
0 


f ‘[adt + jdy + gdu]dt = 0. 
0 


af d 
8 aa $x |dt = 0. 
{a+ 85 y+ g| 0 


Integrating by parts, 


5 dy| — aero Ou: 8 
[+e] - [God - = 
41 /d’x 2y 
= 7 \0 “48 = 


but since 6x and dy are wholly arbitrary, this is impossible unless 
the coefficients of dz and dy in the integrand are both zero. 
ax 
Hence aaa 0, 
oY 
ae 


(6) By the principle of stationary action, we have 


1 eee 
3) ? YP \ ce = 
if m (a + Yaar 0, 


" 


whence ie (Ge + oT art [5 <3@+/)dr=0 early 
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But 7 Et HP) — mgt =h. 
MY “2 
Hence sETP) — go0, 
and s@+ y= sete See 
Bh awene 
edit aha ee rey dy, Gairiee cy anal a 
a = 8e + YY ta w+ os, oy OFG) 7% 


; ; dzx_d dx d 
ee eee 4 : 
since ot Sheet dx aia ov; (vy. App. B; §-6,(1)) 





ea ta d 
ay, S adhe S/N egy 
Eells iia ila a a ee a 


and (1) becomes 


ik (Ca +A qbars [ [eget 95,80 |e 
mie e +P) 5 ear+ [~ gb dr = 0, 


or AL [edartadogors ‘gb = dr = 0. 


Integrating by parts, we get 


| #8e-+ sey] — [7 [epee - dy — y= sr |dr=0 


But dz and dy are zero at the beginning and at the end of 
the actual and the hypothetical paths, so that 


CA 4%, st aH 
i) [@—g) dx + yy] — dr = 9, 
and as 62 and dy are wholly arbitrary, this necessitates that 


og, 


y=, 
as in Art. 47, (a). 
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48. Although we have proved merely that in a system under 
conservative forces the action satisfies the necessary condition 
of a minimum, namely, 64=0, it may be shown by an elabo- 
rate investigation that in most cases it actually is a minimum, 
and that the name “least action,” usually associated with the 
principle, is justified. 

A very pretty corollary of the principle comes from its applica- 
tion to a system moving under no forces or under constraining 
forces that do no work. In either of these cases the potential 
energy V is zero, and consequently the kinetic energy is constant, 
that is, 7= h. , " 

A= f 2 Tat = 2h dt =2h(t,—t,) 
and the action is proportional to the time of transit. Hence 
the actual motion is along the course which occupies the least 
possible time. 

For instance, if a particle is moving under no forces, the 
energy a v’ is constant, the velocity of the particle is constant, 


and as it moves from start to finish in the least possible time, 
it must move from start to finish by the shortest possible path, 
that is, by a straight line. 

If instead of moving freely the particle is constrained to 
move on a given surface, the same argument proves that it 
must trace a geodetic on the surface. 


49. Varying Action.* The action, 4, between two configura- 
tions of a system under conservative forces is theoretically 
expressible in terms of the initial and final codrdinates and the 
total energy (v. Art. 46, (3)), and, when so expressed, Hamilton 
called it the characteristic function. 

In like manner, if 


s=f Lat= (rv de= f ‘(h— 2V) dt, 


* For a detailed account of Hamilton’s method, see Routh, Advanced Rigid 
Dynamics, chap. x, or Webster, Dynamics, § 41. 
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S may be expressed in terms of the initial and final coordinates 
of the system, the total energy, and the time of transit. When 
so expressed, Hamilton called it the principal function. By 
considering the variation produced in either of these functions 
by varying the final configuration of the system, Hamilton 
showed that from either of these functions the integrals of 
the differential equations of motion for the moving system 
could be obtained, and he discovered a partial differential 
equation of the first order that S would satisfy and one that 4 
would satisfy, and so reduced the problem of solving the equa- 
tions of motion for any conservative system to the solution of 
a partial differential equation of the first order. 

It must be confessed, however, that in most cases the advan- 
tage thus gained is theoretical rather than practical, as the 
solving of the equation for S or for A is apt to be at least as 
difficult as the direct solving of the equations of motion. 


CHAPTER V 
APPLICATION TO PHYSICS 


50. Concealed Bodies. In many problems in mechanics the 
configuration of the system is completely known, so that a set 
of codrdinates can be chosen that will fix the configuration 
completely at any time; and the forces are given, so that if 
the system is conservative the potential energy can be found. 
In that case the Lagrangian function Z or the Hamiltonian 
function H can be formed, and then the problem of the motion 
of the system can be solved completely by forming and solving 
the equations of motion. 

In most problems in physics, however, and in some problems 
in mechanics the state of things is altogether different. It may 
be impossible to know the configuration or the forces in their 
entirety, so that the choice of a complete set of codrdinates or 
the accurate forming of the potential function is beyond our 
powers, while it is possible to observe, to measure, and partly to 
control the phenomena exhibited by the moving system which 
we are studying. If from results which have been observed or 
have been deduced from experiment we are able to set up in- 
directly the Lagrangian function, or the Hamiltonian function, 
or the modified Lagrangian function, we can then form our dif- 
ferential equations and use them with confidence and profit. 


51. Take, for instance, the motion considered in the problem 
of the two equal particles and the table with a hole in it (v. 
Art. 8, (@)), and suppose the investigator placed beneath the 
table and provided with the tools of his trade, but unable to 
see what is going on above the surface of the table. He sees 


the hanging particle and is able to determine its mass, its 
98 
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velocity, and its acceleration ; to fix its position of equilibrium ; 
to measure its motion under various conditions; to apply addi- 
tional forces, finite or impulsive, and to note their effects. 

His system has apparently one degree of freedom. It pos- 
sesses an apparent mass m and is certainly acted on by gravity 
with a downward force mg. He determines its position of 
equilibrium, and taking as his codrdinate its distance 2 below 
that position, he painfully and laboriously finds that 2, the accel- 
eration of the hanging mass, is equal to 


dh a 


He is now ready to call into play his knowledge of mechanics. 
If 7 is the kinetic energy and Z is the Lagrangian function for 
the system which he sees, 














T= Te, 
and L=T-V=t#-¥. 
Sj 5 (3) =z 
a _ di\ée) ~ bn’ 
Poaeiee ovVem Ls a’ 
Tip Gp DN" GG Bayt |” 
oVie m ae 
eal Pa ak Sa 
Therefore ame Fe ae g 
Pai | Me Sama, ee |) = 
Se ae | see | we 


The motion, then, can be accounted for as due to the 
downward force of gravity combined with a second vertical 


; : my a : 
force having the potential energy oF Ss +x}, that is, a 


force vertically upward having the intensity 7 ee +1] ; 
and of course this force must be the pull of the string and may 
be due to the action of some concealed set of springs. 
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On the other hand, the moving system may contain some 
concealed. body or bodies in motion, and that this is the fact is 
strongly suggested if a downward impulsive force P is applied 
to the hanging particle when at rest in its position of equilib- 
rium. For such a force is found to impart an instantaneous 


> 





velocity z= es » just half what we should get if the hanging 


am 
particle were the only body in the system, and just what we 
should have if there were a second body of mass m above the 
table, connected with the hanging particle by a stretched string 
of fixed length. 

Obviously this concealed body is ctgnorable, since we have 
already found a function from which we can obtain the differ- 
ential equation for x, namely, the function we have just called 
the Lagrangian function Z. It is 


M 2 mY a ome 
57 Sony a 





and contains the single codrdinate 2. 

But if we have ignored a concealed moving body form- 
ing part of cur system, this expression is not the Lagrangian 
function Z, but is equal or proportional to ®&, the Lagrangian 
function modified for the codrdinate or codrdinates of the con- 
cealed body; and in that case some or all of the terms which 
we have regarded as representing the potential energy of the 
system may be due to the kinetic energy of the concealed body 
(ve Art. 41). 

Of course the complete system may have two or more degrees 
of freedom. Let us see what we can do with two. Take zx 
and a second codrdinate @ and remember that as @ is ignor- 
able it must be a cyclic coordinate and must not enter into 
the potential energy. 

Let us now form the Lagrangian function and modify it for 6. 

We have 7= 42 + Bi + C6”, where A, B, and C are functions 
of 2. 
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oT 
— Bi +2 C6. 





ios 


6 


ay ye \ 2 
T= Av? + oo (gt — Be?) + Cine me : 


To modify for 8 we must subtract 6p, We get 


— Bay’ es Po — Bpot 
4C 2C 





My = Ad’ + 7, Pet — Ba’) + (Pe 


pe prise 2 
= A — ——— Y = g . 
malta 2 5 gree 4C 
Suppose that no external force acts on the concealed particle, 
so that the potential energy of the system is — mgz. 
Then, if the Lagrangian function modified for @ is ®, 


B B 
@ =| 4. |a* 
ag *3¢ 





pot — i + nga. 


Since on our zgnoration hypothesis 6 does not enter the 
potential energy, the momentum ae K, a constant, and 





‘pe Kk? 
® = 4 == leet ai 5 ae Gi nae mMgx. 


But we know that @ is equal or proportional to 


M .. mg a ML 


Deo (geen 92 


the function which on our hypothesis of no concealed bodies 





v 


we called L. 





We see that if B= 0, if d= m, and ee = A=5 ee, 
whence ae : eee 2 
®=2L 
Then we have T= me + ss K*(a wey 62, 
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where K is the momentum corresponding to the codrdinate 6 
and may be any constant. 
If we take for K? the value m’ga’, 


P= FP +e +(a—2)F], 


5 lB +(a— 2/0] 


is obviously the kinetic energy of a mass m moving in a plane 
and having a—-z and @ as polar codrdinates. The Lagrangian 
function is equal to 


5 [22+ (a — v)°6"] + mgn, 


and the x Lagrangian equation is 
2 mi = — m(a — 2) 6 + mg. (1) 


If the angular velocity of the concealed mass when the hang- 
ing particle is at rest in its position of equilibrium is @,, since 


in that case x= 0 and % = 0, equation (1) gives us 6, = ve The 
a 


observed motion of the hanging particle is then accounted for 
completely by the hypothesis that it is attached by a string 
to an equal particle revolving on the table and describing a 


circle of radius a@ about the hole in the table with angular 
if 


velocity A when the hanging particle is at rest in its posi- 
a 


tion of equilibrium. We see that on this hypothesis the term 
3 
eas || which on the hypothesis that the system 


contained only the hanging particle was an unforeseen part of 
the potential energy, is due to the kinetic energy of the con- 
cealed moving body. 

It may seem that giving & a different value might lead to 
a different hypothesis as to the motion of the concealed body 
that would account for the motion of the hanging particle. 
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Such, however, is not the case. We have 





eek (a ay)" 
roe toa ae | (2) 
VG 
Let = sey) 
® mv gas 
Then P= FP +e + (a2) $i). 


“ [27+ (a — x)? ¢’] is, as above, the kinetic energy of a body 


of mass m, with polar codrdinates a — z and ¢; and the concealed 
motion is precisely as before. In using the form (2) we have 
merely used as our second generalized codrdinate @, a perfectly 
suitable parameter but one less simple than the polar angle. 


52. Problems in Physics. In physical problems there may be 
present electrical and magnetic phenomena and concealed 
molecular motions, as well as the visible motions of the mate- 
rial parts of the system. In such cases, to fix the configuration 
of the system even so far as it is capable of being directly 
observed we must employ not only geometrical codrdinates 
required to fix the positions of its material constituents, but also 
parameters that will fix its electrical or magnetic state; and as 
we are rarely sure of the absence of concealed molecular 
motions, we must often allow for the probability that the func- 
tion we are trying to form by the aid of observation and 
experiment and on which we are to base our Lagrangian equa. 
tions of motion may be the Lagrangian function modified 
for the ignored codrdinates corresponding to the concealed 
molecular motions. 


53. Suppose, for instance, that we have two similar, parallel, 
straight, conducting wires, through which electric currents due 
to applied electromotive forces are flowing. It is found ex- 
perimentally that the wires attract each other if the currents 
have the same direction, and repel each other if they have 
opposite directions, and that reversing the currents without 
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altering the strength of the applied electromotive forces does 
not affect the observed attraction or repulsion. It is known 
that an electromotive force drives a current against the resist- 
ance of the conductor, and that the intensity of the current is 
proportional to the electromotive force. Moreover, it is known 
that an electromotive force does not directly cause any motion 
o: the conductor. It is found that, as far as electric currents 
are concerned, the phenomena depend merely on the intensity 
and direction of the currents. 

To fix our configuration we shall take x as the distance 
between the wires and take parameters to fix the intensities of 
the two currents. These parameters might be regarded as 
coordinates or as generalized velocities, but many experiments 
suggest that they are velocities. We shall call them y, and y, 
and define y, as the number of units of electricity that have 
crossed a right section of the first wire since a given epoch. 

As all effects depend upon y, and y, and not on y, and y,, 
y, and y, are cyclic coordinates. 

Let us suppose that there are no concealed motions. Then 
the kinetic energy 7 is a homogeneous quadratic in @, y,, and 
Yor Let G +9 +9 ee 50) ae a0 

T= Ax’ + Ly; + My y, + Ny; + Bay, + Cay, (1) 
where the coefficients are functions of 2. 

Since reversing the directions of the currents, that is, revers- 
ing the signs of y, and y,, does not change the other phenomena, 
it must not affect 7; therefore B and C are zero. 

If the wires are not allowed to move, «=0 and 7 reduces 
to Lys + My,y,+ Ny, which is called the electrokinetic energy 
of the system. Since from considerations of symmetry this can- 
not be altered by interchanging the currents, NV = L. 

Let us now suppose that the first wire is fastened in position 
and that the only external forces are the electromotive forces 
E, and £,, producing the currents in the two wires; the resist- 
ances Ry, and Ry, of the wires, equal, respectively, to #, and 
#, when the currents y, and y, are steady; and an ordinary 


Cuar. V] PARALLEL LINEAR CONDUCTORS 105 


mechanical force F, tending to separate the wires. We now 
have 


and for our 2 Lagrangian equation 





day dA dl, GM. al 
Big 


2 Av + 2 Se Oy PN OF 
“fa d dx da Ki dx Ws dx Une) 


Let us study this equation. First suppose the electromotive 
forces and the impressed force F all zero, so that 9,=9,=%, 
and F=0. Equation (2) reduces to 


Dede eet 0, 
dx 


or SS = ee, 


If, then, a transverse velocity were impressed on the second 

wire, the wire would have an acceleration unless “= 0. But 
aD 

both wires, on our hypothesis, being inert, they can neither 


attract nor repel each other, therefore A is a constant; and as 
T = Ax’ when y,=¥y,=9, A is a positive constant. Therefore 
9 

(2) reduces to ff ya aT 


De elle es Le tad 
AD a Us (8) 


Let us now suppose that y,=0, and F=0. Equation (3) 


becomes 2 At — a Y, = 0, 
da 
or pete ll ie 
2A de?” 


and the second wire is attracted or repelled by the first, 
even when no current is flowing through the second wire. 


ae dL ; 
This is contrary to observation. Therefore [Pe 0, and LZ is 


a constant; and as 7’= Ly? when # and y, are zero, L is a 
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positive constant. Equation (3) now becomes 


24¢—-“ yy, =F, 
: ek 40 
and if F=0, ao fs 


ay Rs aM 
If y, and y, are of the same sign, % and os have the same 
v 


sign. But according to observation the wires attract if the 
currents flow in the same direction, therefore % is negative and 
dM . : 
—— ls negative. 
dx 

If x=0, we have a single wire carrying a current of inten- 
sity y¥,+y,- The electrokinetic energy Ly? + My,y, + Ly; be- 
comes L(y,+¥y,)* or Ly +2Ly,y, + Ly;, so that the value of 
M when x=0 is 2Z. Hence, as M is a decreasing function, M 
is less than 2 Z always. As M is easily seen to be zero when 
x is infinite, it must be positive for all values of a. 

We have, then, 

T= Ad’ + Ly? + My,y,+ Ly, (4) 
where 4 and JZ are positive constants and J is a positive decreas- 
ing function of z always less than 2 LZ. ZL is called the coeffictent 
of self-induction of either wire per unit of length, and M the coeffi- 
cient of mutual induction of the pair of wires per unit of length. 

Our Lagrangian equations are 


OE 
2 Ai ——— 9, = Fr, (5) 
CF Oe ma . ; 
aL? Ly, + My,\= £,— Ry, (6) 
d ; : a : 
er [2 Ly, + My,] = E, — Ry,. Gr) 


54. Induced Currents. (a) Suppose that no current is flow- 
ing in either of the wires considered in the last section, and 
that the first wire is suddenly connected with a battery fur- 
nishing an electromotive force H,, and that thereby a current 
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y, 8 impulsively established. Then, by Thomson’s Theorem 
(v. Art. 30), such impulsive velocities must be set up in the 
system as to make the energy have the least possible value 
consistent with the velocity caused by the applied impulse. 
If the current y, set up in the first wire has the intensity ¢, 
T= Az?” +. It? + My, + Ly}; 

and making Z’ a minimum, we have 

2 At = 0, 

Mi+2Ly,=0; 


whence z= 0, and the second wire wili have no initial velocity, 
Mi 
21 
second wire, of intensity proportional to the intensity of the cur- 
rent in the first wire. Since, as we have seen, Jf and Z are both 
positive, y, is negative, and this so-called induced current will be 
opposite in direction to the impressed current 7,. This impulsively 
induced current is soon destroyed by the resistance of the wire. 

(6) Suppose that a steady current 7,, caused by the electro. 
motive force #,, is flowing in the first wire while the second 
wire is inert, and that #, and consequently ¥y» is impulsively 
destroyed by suddenly disconnecting the battery. This amounts 
to impulsively applying to the first wire the additional electro- 
motive force —#,. If the system were initially inert, this, as 
we have just seen in (@), would immediately set up the induced 


, and a current will be set up impulsively in the 





and ¥,=— 


: MM .. : ; é 
eurrent y, = >— in the second wire, and we should have y, = —1, 


21 
I. = aa as the immediate result of our impulsive action. Combine 
this with the initial motion in our actual problem, y, =7, 7, = 9, 
é Mi 
and we get for our actual result 7, = 0, 7, = ia (v. Art. 34). 
So that if our first wire is suddenly disconnected from the 
battery, an induced current whose direction is the same as that 


of the original current is set up in the second wire. It is, how- 
ever, soon destroyed by the resistance of the wire. 
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(c) Suppose that we have a current g, in our fixed wire, caused 
by a battery of electromotive force #,, and no current in our 
second wire, and that the second wire is made to move away 
from the first. Equations (6) and (7) of the preceding section 
give us 


dy 


(41? — M”) ai = 2L(#,—Ry,)—-M(L,— By,) 


~(2Ly, —My,) sm 

(41° a’) “OLB, — Ry,)— M(E, — By.) 
: OAL 
= Clg Nye 


When we are starting to move the second wire, 


H = 0, y= 0, HE, — Ry, = 9, 
and we have 


} » ay au 
(41? — mM”) Ti UE eee 


dy _.dM 
41? —M?) -2=— 219g ¢-—. 
( sary I” Te 


As we have seen, L, M, 4 L’ — M’, are positive and a is nega- 
x 


tive. Hence the current y, will decrease in intensity, and a 
current y, having the same Ee as y, will be induced in 
the moving wire. 

The phenomena of induced currents which we have just 
inferred from our Lagrangian equations are entirely confirmed 
by observation and experiment. 
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SYLLABUS. DYNAMICS OF A RIGID BODY 


1. D’Alembert’s Principle. In a moving system of particles the 
resultant of all the forces external and internal that act on any 
particle is called the effective force on that particle. Its rectangular 
components are mz, mij, and mz. 

The science of rigid dynamics is based on D’ Alembert’s prin- 
ciple: In any moving system the actual forces impressed and 
internal, and the effective forces reversed in direction, form a set 
of forces in equilibrium, and if the system is a single rigid body, 
the internal forces are a set separately in equilibrium and may 
be disregarded. 

It follows from this principle that in any moving system the 
actual forces and the effective forces are mechanically equivalenv. 
Hence, 

(a) Smé = BX, 

(b) Sm [y# — xy] = S[yX — 2], 

(c) Sm [#dax + dy + 282] = S [Ne + Voy + Zz]. 

These equations may be put into words as follows: 

(2) The sum of those components which have a given direction 
is the same for the effective forces and for the actual forces. 

(6) The sum of the moments about any fixed line is the same for 
the effective forces and for the actual forces. 

(c) In any displacement of the system, actual or hypothetical, the 
work done by the effective forces is equal to the work done by the 
actual forces. 

Equations (a) and (?) are called differential equations of motion 
for the system. 


2. p, = Smv, = mz and is the linear momentum of the system in 
the X direction; h, = Xm[yv,—xv,] = Xm[yt — ay] and is the 
moment of momentum about the axis of Z. 
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Equations (a) and (6) of § 1 may be written, respectively, 





tbe = 3X,- and. 


and § 1, (a), and § 1, (6), may now be stated as follows : 

(2) In a moving system the rate of change of the linear momentum 
in any given direction is equal to the sum of those components of 
the actual forces which have the direction in question. 

(b) In a moving system the rate of change of the moment of 
momentum about any line fixed in space is equal to the sum of the 
moments of the actual forces about that line. 


TL aa 
oe S(ysX —aY |, 








Pa nx 
3. Center of Gravity. z= 2m A 
M 
Hence 0, = aie = - a 


or, the linear momentum in the X direction is what it would be if 
the whole system were concentrated.at its center of gravity. 


a ij 
b, =u |g G24 a4) 4 3mly ee ot, 


or, the moment of momentum about the axis of Z is what it would be if 
the whole mass were concentrated at the center of gravity plus what 
it would be if the center of gravity were at rest at the origin and the 
actual motion were what the relative motion about the moving center 
of gravity really is. 

4. The motion of the center of gravity of a moving system is the 
same as if all the mass were concentrated there and all the actual 
forces, unchanged in direction and magnitude, were applied there. 

The motion about the center of gravity is the same as if the center 
of gravity were fixed in space and the actual forces were unchanged 
in magnitude, direction, and point of application. 


5. If the system is a rigid body containing a fixed axis, 
h, = Mk" = M(W? a k?) w, 


where Mk? = M(h? + k°) and is the moment of inertia about the axis, 
and where w is the angular velocity of the body. 


Equation (0) of § 1 becomes une = N, where N is the sum 


of the moments of the impressed forces about the fixed axis. 
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6. If the system is a rigid body containing a fixed point and o,, 
w,, ,, are its angular velocities about three axes fixed in space and 
passing through the fixed point, 


h, = Cw, — Ew, — Dw,, 
where C is the moment of inertia about the axis of Z, and D and FE 
are the products of inertia about the axes of X and Y, respectively ; 
that is, 
C= ain (a 4-4"), D=Smyea, E = Smzex. 
Equation (/)) of § 1 becomes 


dh dw dw dw 
eet) 2— fF — — p—! Je 18, re 
dt dt dt gp a) SO He, 





— Do,w, — F (02 — o7)=N, 


where NV is the sum of the moments of the impressed forces about 
the axis of Z. 


7. Euler’s Equations. If the system is a rigid body containing 
4 fixed point and o,, w,, ,, are its angular velocities about the 
principal axes of inertia through that point (a set of axes fixed in 
the body and moving with it), equation (/) of § 1 becomes 


8. Euler’s Angles. Euler’s angles 0, w, ¢, are codrdinates of a mov- 
ing system of rectangular axes A, B, C, referred to a fixed system 
X, Y, Z, having the same 
origin O. @ is the colati- 
tude and wy the longitude 
of the moving axis of C 
in the fixed system (re- 
garded as a spherical sys- 
tem with the fixed axis of 
Z as the polar axis), and 
@ is the angle made by 
the moving CA-plane with 
the plane through the fixed 
axis of Z and the moving 
axis of C. 


Z 
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We have o, = 6 sin d — sin 6 cos 4, 
w, = 6 cos ? + ysin G sin ¢, 
o, = y cos 6 + db; 
and w, =— Osing + ¢ sin 6 cos y, 
w, = 6 cosy + ¢ sin 6 sin y, 
w, = $ cos 6 + Wp. 
9. 4 mv?, the kinetic energy of a particle, becomes 
a (a? + 7 + 2] in rectangular codrdinates, 
Ue me ° 26 : 
3 [7 + r°p?] in polar codrdinates, 


es 4 77 (@ + sin?6¢”) ] in spherical codrdinates. 


m ia ; 
10. » 9 v’, the kinetic energy of a moving system, becomes 


™m 
We) 2 PAN _ ard i . 
y 2 (a? + 7 + 2) in rectangular codrdinates ; 


A : eo 5 : 
ko? if a rigid body contains a fixed axis ; 





M[ (dz? aay) 7 Sa 
9 ( =) +(2) + 9 k?w? if the body is free and the motion 


is two-dimensional ; 


$[4o,? + Bw} + Cu,’ — 2 Dw,w, — 2 Fw,w, — 2 Fw,w,] if the 
body is rigid and contains a fixed point ; 

$[Aow/? + Bu,’ + Cw,"] if the body is rigid and the axes are 
the principal axes for the fixed point. 


11. Impulsive Forces. In a system acted on by impulsive forces, 
the resultant of all the impulsive forces external and internal that 
act on any particle is called the effective impulsive force on that 
particle. Its rectangular components are m(#,—4,), m(Y¥, — Y%)s 
(mz, — 2,). D’Alembert’s principle holds for impulsive forces. It 


APPENDIX A 113 


follows that in any system acted on by impulsive forces, the actual 
forces and the effective forces are mechanically equivalent. Hence, 

(a) 2m[a#, — &,]= =X, 

(0) Xm Lye, a xy, oa YX te Yo | = sya — x¥], 

(c) 2m[ (a, — &,) da +(¥, — %) dy + (4, — 4,) 82] 

= &[ Xda + Voy + Z6z]. 

These equations may be put into words as follows : 

(a) The sum of the components which have a given direction 
is the same for the effective impulsive forces and for the actual 
impulsive forces. 

(6) The sum of the moments about a given line is the same for 
the effective impulsive forces and for the actual impulsive forces. 

(c) In any displacement of the system, actual or hypothetical, 
the sum of the virtual moments of the effective impulsive forces is 
equal to the sum of the virtual moments of the actual impulsive 
forces. 

Equations (a) and (/) are the equations for the initial motion 
under impulsive forces. (a) and (4) may be restated as follows: 

In a system acted on by impulsive forces, the total change in the 
linear momentum in any given direction is equal to the sum of those 
components of the actual impulsive forces which have the direction 
in question. 

In a system acted on by impulsive forces, the total change in the 
moment of momentum about any fixed line is equal to the sum of 
the moments of the actual impulsive forces about that line. 

Section 4 holds unaltered for impulsive forces. 


APPENDIX B 
THE CALCULUS OF VARIATIONS 


1. The calculus of variations owed its origin to the attempt to 
solve a very interesting class of problems in maxima and minima 
in which it is required to find the form of a function such that the 
definite integral of an expression involving that function and its 
derivatives shall be a maximum cr a minimum. 

Take a simple case: If y= f(x), let it be required to deter- 

x 
mine the form of the function f, so that si ‘ol ¥, i= ate shall be 
4% maximum or a minimum. I 

Let f(x) and F(x) be two 

possible forms of the function. 


dy 
dx 


Consider their graphs y = f(x) 
and ¥ = F(a). 

If n(x) is F(x) — f(x), n(x) 
can be regarded as the inere- 
ment given to y by changing 
the form of the function from 
f(x) to F(x), the value of the 
independent variable x being 
held fast. 

This increment of 7, 7 (x), is called the variation of y and is written 





Co x v, 


d7; it is a function of 2, and usually a wholly arbitrary function of a, 


Vane i dy 
The corresponding increment in v', where 7' = —-> can be shown to be 
wv 


dy 
n'(v), and is the variation of y',and is written 3, or 8. =n Oby. iously, 
Js 


ad 


$y! 
Ze dix 


Sy. (1) 


Tf an infinitesimal sr dy is given to y, it is proved in the 
differential calculus that - OL differs by an infinitesimal of 


higher order than 8y mee a increment produced in $(y). This 
114 
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approximate increment is called the variation of $(y), so that 


26 (y) =F by) (2) 
Similarly, Baan v) = Fo by + FS ay (3) 
or since a is not varied, 
sate ye as by + 7 Sy). Lae 
As Ib) = 7, 6) 
and db(y, y')= - dy + a dy, 


we can calculate variations by the familiar formulas and processes 
used in calculating differentials. 

2. Let @ be an independent parameter. Then y = f(r) + ay(2), 
or y=f(x)+ aby, is any one of a family of curves including 
y = f(x) (corresponding to a= 0) and y= f(a) + 9(«) (corresponding 
toa=1). If x, and z, are fixed values, and if 


1(@) =f $ (4, y + aby, y' + aby') dz, 
*o 
I(a) is a function of the parameter @ only. A necessary condition 
that /(a), a function of a single variable a, should be a maximum 


o2 a minimum when a = 0 is I'(a@) = (0) when a = 0. 


zy ) 
I'(a) =|) A POEs Y + ay, y' + aby') da 


Te ang (ah a er ] , 
I'(a@) =i eo by + By! by' | dx 
when a= 0. 
That is, 1'(0) = i “8p de. (v. $1, (4)) 


A necessary condition, then, that f b (2, y, y')d« should be a 


Ly ey 
maximum or a minimum when y = f(z) is iy dp (a, y, y')dx = 0. 
ty 
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J eae is taken as the definition of the variation of lh pdx, and 


our necessary condition is usually written 


of 6 (@, GY iti), 


How this condition helps toward the determination of the form 
of f(#) can be seen from an example. 


3. Let it be required to find the form of the shortest curve 
y = f(x) joing two given points («,, y,) and (a, y,). Here, since 


= V14 y"dz, 
I =f ‘VI y®da, 


and / is to be made a minimum. 


ee =I “3VIidy"-d 
mle yoy 
Vi+/° +y” 
f tee y' ° by 
by 9 Vi+y? a Y) ‘ 
tae % 
= il Vityt by. 


Integrating by parts, this last reduces to 


y! es 7rd y' 
SSS oy) — —— Ss oy) dx. 
Ea at y? a i dx V/1 a y” UI 


y" are 
Vipg (ers 


since, as the ends of the path are given, ’y = 0 when «=a, and 
when x=2,. Then J = 0 if 














id Yy 
J. EV oe 


but since our dy (that is, »(#)), is a function which is wholly arbitrary, 
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d ! 
the other factor, — —_ must be equal to zero if the integral 
Y 





: rn 
is to vanish. eee 4 
set 
This gives us le tT 
V1i+y" 
yf 0 
=ca+d; 


and the required curye is a straight line. 


4. In our more general problem it may be shown in like manner that 


of 6 (®, y, y')dx = 0 


leads to a differential equation between y and x and so determines y¥ 
as a function of 2. 

Of course 6/ = 0 is not a sufficient condition for the existence of 
either a maximum or a minimum and does not enable us to diserim- 
inate between maxima and minima, but like the necessary condition 


if 


sams 0 for a maximum or minimum value in a function of a single 
dx 


variable, it often is enough to lead us to the solution of the problem. 
5. Let us now generalize a little. Let «, y, z,---, be functions of 


: : dx dy dz 
an independent variable 7, and let a’ = = 25 Up Sp 
dy dy dr 


Suppose we have a function $(7, 2, y, 2,--+,2', y', ',+++). 
By pre the forms of the functions ae pee r fast, let 
XL, Y, 2++ , be given the increments €(7), y(7), €(7),-- 
AD hod becomes x + &(7r), y becomes 7 + 7(7), # becomes + €(7), ++ +5 
x' becomes a! + €'(7), y' becomes y' + 7'(1), 2’ becomes z'+ €'(r), +++ 
E(r), E(r), ave the variations of x and «' and are written dx 
and $x’. Obviously, d 
Op —— On; 
dy 
The increment produced in ¢ when infinitesimal increments 4.7, 
dy,---, dx’, dy',---, are given to the dependent variables and to their 
derivatives with respect to 7 is known to differ from 
£8 sn + Day 4 Fas 28 bot Bay! + 
Ox dy Ox: 
by terms of higher order than the variations eee This approxi- 
imate increment is called the variation of ¢ and is written 3¢. It is 
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found in any case precisely as dd, the complete differential of ¢, 


is found. 


That f do(7, 2, y,+++, x', y',-++)dr =0 is a necessary condition 


that f~ P(", L, Y,+++,2', y',++-) dr should be a maximum or a 
minimum can be established by the reasoning used in the case of 
Ae yy yan. “Lhe integral f NE dr is called the variation 
ry a re 
of ff odr, so that af pdr =f 8 dr. 


6. It should be noted that our important formulas 


and 5f ‘gar = “8 dr 


hold only when r is the independent variable which is held fast 
when the forms of the functions are varied; that is, when 67 is 
supposed to be zero. 

If x and y are functions of 7 and we need ait, we get it in- 
directly thus: oe 





ay _y 
dx x! 
Dy EU ee ee ey =| a8 dy d |, 
OR wa) Sa a ara ~~ gt! ‘ine dx Fie 
dy — d dy d 
pe tae ° dx dx oi dx, dz or: > 


If we need 6 iI gdx, we must change our variable of integration 


to 7. oe 
eae = | oar 
Bf gde=s fo ore 
Bf eae = [39 (+5 ale (2) 
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